STEADY  STATE  STRESS  IN  A  COATED 
INFINITE  HALF-SPACE  SUBJECTED  TO  A 
MOVING  LOAD 


THESIS 

Jason  M.  Cruthirds,  1st  Lieutenant,  USAF 
AFIT/GCS/ENC/05-0 1 

DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 

Wright-Patterson  Air  Force  Base,  Ohio 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  of  the  United  States  Air  Force,  Department  of  Defense,  or  the  United  States 
Government. 


AFIT/GCS/ENC/05-0 1 


STEADY  STATE  STRESS  IN  A  COATED  INFINITE  HALF-SPACE 
SUBJECTED  TO  A  MOVING  LOAD 

THESIS 


Presented  to  the  Faculty 
Department  of  Mathematics  and  Statistics 
Graduate  School  of  Engineering  and  Management 
Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Computer  Systems 


Jason  M.  Cruthirds,  BS 
1st  Lieutenant,  USAF 

March  2005 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


AFIT/GCS/ENC/05-0 1 


STEADY  STATE  STRESS  IN  A  COATED  INFINITE  HALF-SPACE 
SUBJECTED  TO  A  MOVING  LOAD 


Jason  M.  Cruthirds 
1st  Lieutenant,  USAF 


Approved: 


//Signed//  1 1  March  2005 

William  P.  Baker  (Chainnan)  date 

//Signed//  1 1  March  2005 

Anthony  N.  Palazotto  (Member)  date 

//Signed//  1 1  March  2005 


Lawrence  K.  Chilton  (Member) 


date 


AFIT/GCS/ENC/05-0 1 


Abstract 

This  research  investigates  the  use  of  coatings  to  mitigate  the  stress  distribution 
into  an  infinite  half-space.  High  energy  impact  phenomenon  at  velocities  exceeding  the 
speed  of  sound  is  an  important  area  of  research  to  the  Air  Force  Research  Laboratory. 
Holloman  Air  Force  Base’s  High  Speed  Test  Track  sustains  significant  damage  due  to 
this  phenomenon.  In  this  thesis,  the  track  system  and  coating  are  modeled  analytically 
with  equations  of  motion  in  terms  of  linear  displacements.  Coating  thickness  and 
material  properties  of  epoxy  or  polymer  laminates  are  investigated  to  understand  their 
affect  on  stress  distribution  in  the  rail.  An  analytic  solution  is  used  to  verify  the 
numerical  solutions.  It  is  found  that  due  to  the  limitations  in  coating  thickness  of  the 
track  system,  this  property  has  no  significant  affect  on  the  stress  distribution.  However, 
the  shear  modulus  of  the  material  is  found  to  have  a  significant  affect  representing  the 
possible  onset  of  material  failure  through  the  consideration  of  the  combined  stress  field  at 
the  coating-rail  interface. 
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STEADY  STATE  STRESS  IN  A  COATED  INFINITE  HALF-SPACE 
SUBJECTED  TO  A  MOVING  LOAD 

Chapter  1.  Introduction 

1.1  General  Issue 

High  energy  impact  phenomenon  at  velocities  exceeding  the  speed  of  sound  is  an 
important  area  of  interest  to  the  Air  Force  today.  Specifically,  the  Air  Force  Research 
Laboratory’s  Holloman  AFB  High  Speed  Test  Track  (HHSTT)  has  a  special  interest  in 
studying  this  phenomenon. 

The  HHSTT  has  a  rocket  sled  attached  to  a  rail  by  a  slipper.  The  rocket  sled  is  a 
moving  vehicle  used  to  obtain  velocities  of  exceptional  magnitude  in  order  to  test 
propulsion  components  of  interest  to  the  Air  Force.  On  April  29,  2003,  HHSTT’s 
hypersonic  upgrades  allowed  a  Missile  Defense  Agency  payload  to  set  a  world  land 
speed  record  at  Mach  8.5  (2.8  km/sec).  The  goal  of  the  HHSTT  is  to  operate  up  to  Mach 
10  (3  km/sec)  in  air  (1:1-1). 

During  high  speed  impact  tests,  the  rail,  composed  of  1080  steel,  sustains 
significant  damage  from  a  phenomenon  known  as  hypervelocity  gouging.  Since 
repairing  the  damage  of  the  rail  is  such  a  great  expense,  the  HHSTT  and  the  Air  Force 
Office  of  Scientific  Research  seek  methods  to  mitigate  gouging.  One  method  currently 
used  is  to  paint  thin  polymer  coatings  on  the  rail.  However,  very  little  is  known  about  the 
affect  of  these  coatings  on  the  stresses  in  the  steel  rail  that  lead  to  gouging. 

1.2  Research  Objectives 

This  research  investigates  the  use  of  coatings  to  mitigate  stress  wave  propagation 
into  a  half-space  by  high  energy  created  by  a  load  moving  parallel  to  the  surface  at 
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constant  hypersonic  velocities.  Particularly,  it  addresses  the  issues  of  coating  thickness 
and  shear  modulus  and  their  affect  on  the  stress  distribution  within  the  rail.  Further,  it 
examines  variable  material  properties  through  the  thickness  of  the  coating  to  analyze  the 
effectiveness  of  layering  coatings.  This  is  carried  out  through  a  steady  state  solution. 

1.3  Scope 

Although  it  is  of  primary  concern  to  the  HHSTT,  this  study  does  not  focus  on  the 
gouging  phenomenon  directly.  There  are  many  factors  that  cause  gouging.  One  of  the 
primary  factors  of  concern  is  the  distribution  of  stress  in  the  rail.  This  research 
investigates  the  development  of  stress  in  a  coated  half-space  prior  to  the  onset  of 
hypervelocity  gouging  to  determine  coating  properties  that  may  mitigate  gouging  before 
the  phenomenon  occurs. 

There  are  many  ways  to  model  Holloman’s  track  system.  This  research  takes  an 
analytic  approach  to  the  problem.  The  rail  is  modeled  mathematically  as  an  elastic  half¬ 
space  with  a  thin  elastic  layer  bound  to  the  half-space  for  the  coating.  The  rail  is 
considered  infinite  in  length  due  to  its  size  relative  to  the  slipper.  Displacement 
equations  are  used  to  describe  stress  distribution  in  both  the  rail  coating  and  the  half¬ 
space  (rail).  Material  properties  of  the  coating  region  are  varied  to  investigate  the  affect 
these  properties  have  on  the  distribution  of  stress  in  the  half-space. 

In  order  to  represent  a  complete  problem  description,  one  would  have  to  consider 
the  nonlinear  wave  equations  that  permit  the  development  of  shock  waves.  This  is  not 
considered  in  this  research,  but  certain  important  parameters,  thickness  and  stiffness, 
related  to  coatings  may  be  characterized  based  upon  a  steady  state  solution.  To  make  the 
problem  linear  and  mathematically  tractable,  several  additional  assumptions  are  made. 
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First,  the  effect  of  friction  is  introduced  through  the  shear  stress  at  the  top  boundary  of 
the  coating.  The  materials  for  the  rail  and  coating  are  assumed  to  be  isotropic, 
homogeneous,  and  obey  Hooke’s  law.  Also,  the  distribution  of  stress  is  assumed  uniform 
through  the  width  of  the  rail.  Thus,  the  system  is  only  modeled  in  two  dimensions  under 
a  plane  strain  assumption.  In  other  words,  this  study  assumes  zero  displacement  through 
the  width  of  the  rail.  Further,  this  thesis  only  allows  for  small  displacements  to  take 
place  and  ignores  the  high  energy  concerns  which  implies 
y,  forces  =  mass  x  accelration  . 

Figure  1  gives  the  geometry  of  the  system  taking  into  account  the  assumptions 
made  for  the  problem.  Region  I  is  the  thin  coating  and  region  II  is  the  rail.  The  length  of 
the  slipper  is  21  and  the  thickness  of  the  coating  is  h  .  The  segment  of  the  coating  where 
the  load  is  applied  is  from  -/  to  / .  The  arrow  pointing  in  the  negative  y  direction 
signifies  the  direction  of  the  direct  load.  The  arrows  pointing  in  the  positive  x  direction 
signifies  the  direction  of  the  shear  stress  and  the  motion  of  the  shoe  traveling  at  speed  cs . 


Figure  1.  Geometry  of  the  Coated  Rail  System 

1.4  Literature  Review 
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The  transportation  industry  is  the  primary  catalyst  for  research  in  half-space  stress 
analysis.  Many  articles  are  published  motivated  by  the  study  of  highways  and  railways, 
with  particular  interests  in  soils  and  pavements.  Although  the  focus  of  this  research  is 
not  on  soils  and  pavements,  much  of  the  theory,  formulation,  and  techniques  used  to 
describe  and  solve  the  systems  apply  directly  to  this  effort.  This  section  outlines  the 
ideas  that  have  contributed  to  this  research  effort. 

1.4.1  Beams 

Beam-type  structures  are  studied  quite  frequently  with  respect  to  half-spaces  and 
moving  loads.  Lu  Sun  studied  the  dynamic  response  of  a  flexible  beam  resting  on  an 
elastic  foundation  (2).  The  governing  equations  are  written  in  terms  of  small 
displacements  and  Green’s  functions  are  found  using  a  spatial  Fourier  transform.  Beam 
response  is  given  in  terms  of  convolutions  with  the  Green’s  functions  in  accordance  with 
linear  operator  theory.  In  another  article,  Sun  studies  Bernoulli-Euler  beams  resting  on  a 
viscoelastic  foundation.  The  foundation  is  modeled  as  a  Winkler  foundation.  Again  Sun 
finds  Green’s  functions  and  presents  the  analytic  solution  to  the  problem. 

H.P.  Lee  studied  the  response  of  a  beam  with  a  moving  source  using  the  equations 
of  motion  for  an  Euler  beam  of  fixed  length  (4).  Lee  nondimensionalized  his  parameters 
and  used  a  Lagrangian  approach  to  develop  numerical  results.  Like  Sun,  Rao  studied 
dynamic  response  of  Bernoulli-Euler  beams  (5).  Rao  solved  the  dimensionless  equations 
using  a  perturbation  method,  expanding  modal  coordinates  in  a  series  taking  advantage  of 

a  small  parameter  in  his  study.  Gbadeyan  and  Oni  (6)  and  Choros  and  Adams  (7) 
performed  similar  studies  on  beam-type  structures. 
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1.4.2  Coatings 

Soldatenkov  discusses  a  thin  elastic  strip  connected  to  an  elastic  half-space  (8). 

He  studied  the  effects  of  a  disk-shaped  load  and  fonnulated  his  problem  in  terms  of  wear 
and  contact  pressure.  Soldatenkov  uses  the  method  of  successive  approximations  to 
determine  the  wear  in  the  strip.  Aleksandrov  and  Arutiunian  study  two  different  types  of 
coating,  a  thin  fluid  layer  and  a  stiffener  (9).  Their  equations  are  formulated  in  terms  of 
mean  displacement  and  tensile  defonnation  and  the  boundary  conditions  are  given  in 
tenns  of  stress.  A  Fourier  transfonn  is  used  to  find  solutions  and  asymptotic  methods  are 
used  to  detennine  approximate  solutions.  Finally,  de  Barros  and  Luco  study  the  steady 
state  response  of  a  layered  viscoelastic  half-space  to  a  moving  point  source  in  three 
dimensions  (TO).  The  moving  load  is  characterized  by  the  distribution  of  body  forces  per 
unit  volume.  The  equations  studied  are  formulated  in  terms  of  displacements  and 
stresses,  de  Barros  and  Luco  use  Fourier  analysis  to  find  solutions  and  a  Fast  Fourier 
Transform  to  interpret  the  results  in  the  spatial  domain. 

1.4.3  Equations  of  Motion 

One  thing  that  varies  greatly  from  study  to  study  is  the  equations  used  to  solve  the 
problems.  Kennedy  and  Herrmann  (11),  Vigak  (12),  Dieterman  and  Metrikine  (13)  and 
Dal’  (14)  use  Airy  stress  functions  which  are  fonnulated  in  tenns  of  displacement 
potentials.  Some  equations  are  used  for  very  specific  applications.  Hu  and  Hartley 

formulate  their  problem  in  tenns  of  plate  deflection  since  they  studied  the  separation 
between  plates  and  half-spaces  (15). 

Verruijt  and  Cordova  use  the  same  form  of  the  displacement  equations  as  this 
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study  (16).  However,  they  model  the  half-space  as  a  viscoelastic  material  rather  than  a 
pure  elastic  material.  They  formulate  the  boundary  conditions  in  tenns  of  stresses. 
Dimensionless  parameters  are  detennined  to  simplify  the  system  and  a  moving 
coordinate  system  is  introduced  to  ensure  source  of  moving  load  is  always  at  the  origin. 
Fung  solves  a  similar  problem  with  a  point  source  (17:258-269).  His  equations  are 
solved  in  terms  of  displacement  potentials.  Fung’s  problem  is  not  the  same  exact 
problem  as  this  research,  but  his  results  are  very  similar  and  comparable. 

1.5  Document  Overview 

Chapter  2  of  this  document  examines  the  theory  and  equations  governing  the 
mechanics  of  the  system.  The  equations  of  motion  are  introduced.  Stress  and  strain  are 
explained  and  used  to  determine  appropriate  boundary  and  initial  conditions.  Further,  a 
coordinate  transfonnation  is  used  to  produce  a  moving  coordinate  system  and  the 
equations  of  motion  are  scaled  to  create  dimensionless  parameters  and  variables.  Finally, 
the  steady  state  system  is  introduced. 

Chapter  3  discusses  the  methods  for  setting  up  the  system  and  the  techniques  for 
solving  the  problem.  It  starts  with  a  discussion  of  the  solutions  to  the  problem  if  no 
coating  were  used.  The  analytic  solution  is  found  with  no  coating  in  the  Fourier 
transfonn  domain  and  the  transform  is  inverted  in  terms  of  stresses.  Since  the  problem 
with  variable  coating  properties  must  be  solved  numerically,  the  analytic  solution  for  the 
problem  with  no  coating  is  presented  in  the  transform  domain  but  is  primarily  used  to 
validate  the  finite  difference  method.  The  finite  difference  method  is  derived  and 
compared  to  the  analytic  solution  in  the  transfonn  domain. 

The  fourth  chapter  starts  with  an  analysis  of  the  problem  where  the  rail  is  directly 
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subjected  to  the  moving  load.  A  point  source,  a  uniform  distribution  of  applied  force, 
and  a  parabolic  distribution  of  applied  force  are  used  to  study  the  behavior  of  stress  in  the 
rail.  Then  a  comparison  is  made  between  applying  the  load  directly  to  the  rail  and  to  a 
coated  rail.  The  thickness  and  shear  modulus  are  varied  to  determine  their  affect  on  the 
stress  distribution  in  the  rail.  Finally,  the  shear  modulus  of  the  coating  is  varied  as  a 
function  of  thickness  to  examine  the  effects  of  layered  coatings  and  a  vanishing  shear 
modulus  at  the  top  of  the  coating. 

Chapter  5  provides  a  summary  of  this  thesis.  First  the  theory,  methods,  and 
solutions  are  reviewed.  Then  a  synopsis  of  the  behavior  of  stress  in  the  rail  with  and 
without  coating  is  presented  along  with  the  significant  findings  in  this  study.  Finally, 
suggestions  for  future  research  are  given. 
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Chapter  2.  Theory  and  Governing  Equations 


2.1  Overview 

This  chapter  discusses  the  general  theory  and  presents  the  equations  which  model 
the  stress  distribution  within  the  rail  and  the  coating.  First,  the  equations  of  motion  are 
given  which  leads  to  a  discussion  on  stress  and  strain  expressed  in  terms  of  displacement. 
Then,  stress  waves  are  discussed  as  an  important  part  of  understanding  the  dynamics  of 
the  system.  Next,  the  continuity,  initial,  and  boundary  conditions  are  detennined.  The 
equations  of  motion  and  the  conditions  are  transfonned  through  a  coordinate 
transfonnation  and  scaled  to  give  dimensionless  equations  of  motion.  Finally,  the  steady 
state  equations  are  presented  to  simplify  the  problem. 

2.2  Equations  of  Motion 

The  governing  equations  of  motion  will  be  expressed  in  terms  of  local 
displacements.  That  is,  a  particle  at  position  (x,  y)  is  displaced  under  loading  to  a 
position  (x  +  u(x, y,  t ),  y  +  v(x,  y,  t ))  as  illustrated  in  Figure  2. 


y 


(x+u,  y+v) 


i - ►  x 

Figure  2.  Example  of  Particle  Displacement  in  Terms  of  Local  Displacements 


The  functions  u  and  v  are  used  to  construct  the  system  of  differential  equations  that 
describe  the  behavior  of  particles  in  the  coating  (region  I)  and  the  rail  (region  II)  (18:486) 
(See  Figure  1).  Under  the  assumption  of  plane  strain,  the  equations  for  the  coating 
(region  I)  are 


Glg(y)Au1  + 


Gjgjy)  ggi  _  d\ 

1-2l>j  dx  1  dt2 


G,gWAv1+^0^  =  A^L 

1  1  1-2 ux  dy  1  dt 2 


2.2.1 


2.2.2 


and  inside  the  rail  (region  II)  are 


.  G,  de1  d^u-, 

G0  Aw,  + - - - -  =  p„  — t" 

2  2  1-2u2  2  dt2 


G !2  de2  d2v2 

G-,  Av0  + - = - -  =  A — f- 

1-2 u2  8y  dt2 


2.2.3 


2.2.4 
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with 


e,  = 


ft,‘+Al;  *  =  1.2 


<9x  fiy 


2.2.5 


—  oo  <  x  <  oo,  y  <  0,  r  >  o 

where  G  is  the  modulus  of  elasticity  in  shear,/?  is  the  density,  u  is  Poisson’s  ratio  for 
the  material,  and  the  subscript,  k  ,  is  used  to  designate  the  region.  The  symbol  A  is  the 

Laplacian  operator  such  that  An  =  .  The  function  g(y)  is  introduced  to 

account  for  variation  in  the  shear  modulus  through  the  thickness  of  the  coating  and  can 
vary  from  zero  to  one. 

2.3  Stress  and  Strain 

Stress  is  the  force  per  unit  area  as  a  result  of  application  of  a  load  and  can  be 
expressed  in  two  components,  direct  and  shear.  The  direct  stress  at  a  point  is  the 
component  acting  normal  to  any  plane  perpendicular  to  the  reference  axis  while  the  shear 
stress  is  acting  parallel  to  that  plane.  This  document  uses  au  and  rjj  to  denote  direct  and 

shear  stress,  respectively,  where  i  is  the  reference  axis  and  j  is  the  direction  of  the  stress 
component.  Figure  3  shows  the  components  of  stress  due  to  an  applied  force  on  a  two- 
dimensional  reference  cell.  The  arrow  labeled  crvv  is  parallel  to  the  x  axis  and  signifies 

the  direct  stress  perpendicular  to  the  y-z  plane  while  the  arrow  designated  rvv  is  a  shear 


stress  and  acts  parallel  to  the  y-z  plane  in  the  direction  of  the  y  axis. 
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Figure  3.  Stress  Components  in  Two  Dimensions 


The  stress  is  related  to  strain  through  the  constitutive  equation.  The  unit 
elongation  or  direct  strain,  stj ,  is  the  deformation  of  a  body  as  a  result  of  the  direct  stress 


in  the  i  direction.  Likewise,  the  shearing  strain,  ytj ,  is  the  defonnation  as  a  result  of 

shear  stress.  Shear  strain  can  be  thought  of  as  the  angle  of  defonnation.  The  equations  of 
motion  are  used  to  find  displacements.  Strain  can  be  expressed  in  terms  of  derivatives  of 
displacement.  The  following  equations  are  the  components  of  strain  in  a  two- 
dimensional  system  (18:5-7): 


du  dv 

8  =  —  ;  £  =  — 

“  dx  ”  dy 


y  = - r 

/  xy  ^ 


du  dv 
dy  dx 


2.3.1 


2.3.2 


Given  the  components  of  strain,  the  components  of  stress  can  be  derived  using  Hooke’s 
Law  and  are  expressed  as 


f 


cr  =  2  G 


o 


1-2  v 


- e  +  e . 


2.3.3 


11 


f 


V 


2.3.4 


A,  =  2C 


1-2  u6  +  £yy 


(7  =V\(7  +cr  I 

zz  \  xx  yy ) 


T  =  Gy  2.3.5 

xy  /  xy 

with  e  defined  in  2.2.5.  Once  functions  for  displacement  are  found,  these  equations  are 
used  to  detennine  functions  for  stress  and  strain  (18:5-7). 

2.4  Stress  Waves 

When  a  load  is  applied  to  a  medium,  the  material  in  the  medium  acts  to  relieve  the 
pressure  of  the  load.  This  action  by  the  material  to  redistribute  the  pressure  produces  a 
wave  called  a  stress  wave.  The  speed  of  propagation  for  a  stress  wave  is  based  on  the 
material  properties  of  the  loaded  region  (18:485). 

Unlike  air,  an  elastic  medium  can  be  characterized  by  two  types  of  waves, 
distortion  and  dilatation  waves.  Waves  of  distortion,  or  equivoluminal  waves,  produce 
defonnation  assuming  the  volume  expansion  is  zero.  Distortion  waves  consist  of  particle 
shearing  and  rotation  only.  Waves  of  dilatation,  or  irrotational  waves,  produce  volume 
expansion.  The  speed  of  each  wave  type  is  characterized  in  terms  of  material  properties 
without  regard  to  the  amount  of  force  applied  to  the  loaded  region. 

The  stress  wave  speeds  occur  naturally  in  the  equations  of  motion.  For  example, 


in  2.2.3  and  2.2.4,  if  the  equations  are  divided  through  by  G2 ,  the  ratio 


A 

G2 


appears, 


which  is  the  inverse  of  the  square  distortion  speed,  — ,  which  leads  to  the  formula 

c2 
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Now,  the  dilatation  wave  speed,  call  it  cd ,  is  directly  related  to  the  distortion  speed  for 
that  material  in  the  following  manner: 

U  +  2G 

c"l~ 

where  G  is  the  shear  modulus,  A  is  the  Lame’s  constant,  and  p  is  the  density  of  the 
region  of  interest  (18:487).  The  dilatation  speed  is  not  directly  used  in  the  formulation  of 
the  equations  or  the  calculations  of  stress  and  strain  in  this  document,  since  it  can  be 
shown  to  be  proportional  to  the  distortion  speed  (18:490).  However,  the  distortion  wave 
speed  is  inherent  in  the  scaled  equations  of  motion  presented  in  a  later  section. 

2.5  Continuity  Conditions 

At  the  interface  between  the  rail  (region  II)  and  the  coating  (region  I),  several 
continuity  conditions  must  be  satisfied.  First,  the  materials  are  considered  to  be  perfectly 
bonded  and  maintain  that  bond  through  time.  This  implies  we  must  have  continuity  of 
displacement  for  all  x  and  time,  t ,  or 


ux(x,0,t)  =  u2(x,0,t) 

2.5.1 

© 

ii 

■A. 

© 

2.5.2 

The  perfect  bond  assumed  between  the  coating  and  the  rail  also  leads  to  the  continuity  of 

shear  and  direct  stress  yielding 

© 

II 

© 

2.5.3 

(7yySX’0’t)  =  (Tyy2  (X,0,t) 

2.5.4 
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2.6  Initial  and  Boundary  Conditions 


The  state  of  the  system  must  be  defined  at  t  =  0 .  This  study  assumes  quiescent 
initial  conditions  implying  u(x,  y,  0)  =  v(x,  y ,  0)  =  0  and  ut  (x,  y,  0)  =  vt  (x,  y ,  0)  =  0 ,  where 

nt  is  the  first  time  derivative.  As  discussed  in  the  previous  chapter,  the  rail  is  considered 
infinite  in  length  and  is  modeled  as  a  half-space.  This  suggests  a  description  of  stress 
conditions  as  x  and  y  approach  infinity  must  be  stated.  This  study  assumes  stress 

vanishes  at  infinity  which  leads  to  the  conditions:  cr^ ,  axu  — »  0  as  |x|  — >  oo  and 
cr  ,  — ^  0  as  v  — >  -oo . 

yyi  y 

The  impact  load  is  introduced  through  boundary  conditions  applied  at  the  top  of 
the  coating,  y  =  h.  The  load  is  given  in  the  fonn  of  stress  as  a  function  of  position  and 
time: 


Tm  (x,  h,t)  =  f(x  -  cst,  t)H{l  -X- cst ) 

2.6.1 

(7m  (x,  h,  t)  =  -k(x  -  cst,  t)H(l  -  X  -  cst\) 

2.6.2 

—  oo  <  x  <  oo,  t>  0 

where  H(l  -  x  -  c  d)  is  a  Heaviside  function  piecewise  defined  as 


H(l-  x  —  cst ) 


Jo  |x-c/|>/ 

Jl  |x-cst|</ 


This  function  represents  the  load  moving  at  a  velocity,  cs ,  in  the  direction  of  the  positive 


x  axis  (see  Figure  1).  It  acts  as  a  switch  to  turn  the  loading  off  outside  the  region  of  the 
slipper  of  length  21 . 
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2.7  Coordinate  Transformation  and  Scaling 

2. 7. 1  Moving  Coordinate  System 

This  study  is  primarily  concerned  with  stress  near  or  in  the  neighborhood  of  the 
moving  load.  For  this  reason,  a  moving  coordinate  system  is  introduced  so  the  source  is 
always  centered  about  the  origin.  Verruijt  and  Cordova  (16),  Kennedy  and  Herrmann 
(11),  and  de  Barros  and  Luco  (10)  use  the  same  technique  to  ensure  their  moving  loads 
stay  near  the  origin.  To  ensure  the  proper  frame  of  reference,  the  coordinate 
transfonnation  E,  =  x  -  c,  t  is  introduced.  Before  the  transformed  equations  are  presented, 
the  equations’  dimensionless  parameters  are  introduced  and  the  equations  are  scaled. 

2. 7.2  Motivation  for  Scaling 

Many  mathematical  problems  inherently  contain  very  large  or  small  parameters. 
To  take  advantage  of  the  size  of  parameters  and  to  understand  the  contribution  of  each 
tenn  relative  to  the  size  and  strength  of  the  parameters  in  the  system  of  equations,  it  is 
necessary  to  introduce  dimensionless  parameters.  In  nondimensionalizing  equations, 
reference  values  or  “scales”  are  selected  that  are  natural  or  “intrinsic”  to  the  system 
(19:21 1).  This  system  has  three  fundamental  dimensions.  They  are  mass  (M),  length 
(L),  and  time  (T).  All  of  the  dimensioned  quantities  for  the  equations  of  motion  can  be 
written  in  tenns  of  these  fundamental  dimensions.  Table  1  shows  the  parameters  and  the 
values  of  interest  for  this  problem  along  with  their  fundamental  dimensions.  It  is 
necessary  to  know  these  values  to  understand  why  certain  scales  are  chosen. 


15 


Parameter  Name 

Symbol 

Range 

Units 

Fundamental 

Dimension 

Density 

A 

Pi 

1000-1500 

7872 

Kg/m3 

m'l'V 

Shear  Modulus 

g2 

1.3x109-2.4x109 

8xl010 

Kg/ms2 

m‘l'‘t'2 

Shoe  Speed 

Cs 

2000-3000 

m/s 

mVT'1 

Distortion  Speed 

cx 

C2 

1000-1300 

3188 

m/s 

Coating 

Thickness 

h 

1.5x1  O'4- 1.8x1 0'3 

m 

m^'t0 

Shoe  Length 

21 

0.2 

m 

m^'t0 

Poisson’s  Ratio 

A 

°2 

0.3-0.45 

0.29 

dimensionles 

s 

m°l°t° 

Table  1.  Dimensioned  Parameter  Values  with  Range,  Units  and  Fundamental  Dimension 


2.7.3  Buckingham  Pi/Dimension  Analysis 

In  this  section,  theory  is  presented  to  show  that  finding  the  solutions  to 
dimensionless  equations  of  motion  is  equivalent  to  finding  the  solution  to  the  original 
system.  Then,  using  this  theory,  the  dimensionless  parameters  are  found  that  are  natural 
to  the  system  of  differential  equations,  which  leads  to  the  scaled  equations  in  the 
following  section. 

Let  /  (ql,q2,...,qm )  =  0  be  a  unit  free  physical  law,  meaning  that  it  is  independent 
of  the  particular  units,  like  British  or  metric,  and  {<?,  }”,  be  dimensioned  quantities.  Also, 

n 

let  {Z.}”  be  fundamental  dimensions  with  [q:]  =  ]~ [  Lf  for  all  i  =  \,2,...,m,  and D  be 

1  j= i 

an  nxm  dimension  matrix  with  entries  dji .  If  D  has  rank  r  then  there  are  k  =  m-r 
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independent  dimensionless  quantities,  nx,n2,...,nk  ,  which  are  formed  from  jf/,.}™  .  The 

Buckingham  Pi  Theorem  states  that  the  law  / (qi,q2,...,qm)  =  0  is  equivalent  to 

F{nx,n2,...,n^)  =  0  expressed  in  terms  of  the  dimensionless  quantities  (19: 195-221). 

Through  dimension  analysis  the  natural  scales  can  be  found  by  detennining  the 
dimensionless  quantities  for  the  system. 

For  this  problem,  the  dimensioned  quantities  are  {<jr.}7  =  {pl,p2,Gl,G2,cs,h,l} . 
Let  n  =  [p,  ]“'  [ p2  p  [Gjf  [G2  ]“‘  [cj“5  [/z]“6  [/]“7  =  M if2  T/l' ,  where  n  is  any  quantity 
made  up  of  multiple  powers  of  [qt]  and  [  ]  mean  the  dimension  of  the  listed  quantity. 
Then  define  (3  =  Da ,  where  a  and  (3  are  vectors  of  powers  from  the  definition  of  n . 
With  this  definition,  a  must  be  in  the  null  space  of  D  for  n  to  be  dimensionless.  The 
dimension  matrix  for  this  problem  is 


\P\  Pi  ^2  h  /| 


1 

1 

1 

1 

0 

0 

o' 

M 

D  = 

-3 

-3 

-1 

-1 

1 

1 

1 

< 

L 

0 

0 

-2 

-2 

-1 

0 

0 

r 

When  the  distortion  speeds  are  removed,  the  columns  of  D  correspond  to  the 
dimensioned  quantities  that  occur  in  the  differential  equations  and  boundary  conditions  in 
the  same  order  as  Table  2.1.  The  rows  of  D  correspond  to  the  fundamental  dimensions, 
mass,  length,  and  time,  respectively.  The  vectors  in  the  braces  are  not  part  of  the 
matrix  D ,  but  are  references  used  to  relate  the  entries  in  D  to  the  actual  parameters. 
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Now  rank(D)  -  3  which  implies  there  are  k  =  7-3  =  4  independent  dimensionless 


quantities,  tt,  ,  tt2  ,  ,  nA . 

Now  to  find  the  dimensionless  quantities,  four  linearly  independent  vectors  are 
chosen  from  the  null  space  of  D  .  Table  2  shows  the  vectors  chosen  with  the 
dimensionless  quantity  produced. 


Null  Space  Vector 

Dimensionless 

Quantity 

(0, 0,0,0, 0,1, -1) 

h/1 

^2 

*3 

(0,0, -1,1, 0,0,0) 

G;/G, 

*4 

(-1,1, 0,0, 0,0,0) 

Pi/  A 

1 

fable  2.  Independent  Dimensionless  Quantities 

Since  -  4 ,  all  other  dimensionless  quantities  can  be  written  as  a  product  of 

powers  of  7i^7i1,7i^,nA. 

2. 7. 4  Scaling 

First  an  appropriate  length  scale  is  determined.  With  the  presence  of  the  thin 
coating,  it  is  necessary  to  scale  all  parameters  in  the  y  direction  by  the  thickness  of  the 
coating,  h  .  However,  the  x  direction  is  scaled  by  half  the  size  of  the  slipper,  / .  The  rail 
is  modeled  as  a  half-space  and  has  no  natural  length  scales.  It  will  be  scaled  in  the  same 
manner  as  the  coating  region.  Scaling  the  lengths  in  this  manner  leads  to  the  small 
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rl  3 

dimensionless  parameter,  8  =  nx=  —  ,  which  is  on  the  order  of  10  and  can  give  insight 

into  how  terms  are  contributing  to  the  equations  of  motion. 

The  fundamental  dimension  of  mass  occurs  naturally  in  the  density  and  shear 
modulus  of  the  materials.  However,  from  2.4.1,  it  is  seen  that  mass  and  density  are 
related  through  the  distortion  speed.  Further,  by  dividing  the  equations  of  motion,  2.2.1, 
2.2.2,  2.2.3,  and  2.2.4  through  by  their  respective  shear  modulus,  G, ,  the  inverse  of  the 

naturally  appears.  For  these  reasons,  mass  is  scaled 

through  the  ratio  of  density  to  shear  modulus  introducing  the  speed  scale,  c2 ,  and  given 
the  material  properties  of  both  regions,  a  dimensionless  speed  is  produced  in  the  ratio  of 

shoe  speed  to  distortion  speed,  c  =  =  —  . 

c2 

There  is  another  speed  which  occurs  naturally  and  must  be  considered;  the  slipper 

speed  or  speed  of  the  load,  cs .  The  slipper  speed  is  more  directly  introduced  into  the 

equations  of  motion  through  the  coordinate  transformation,  g  =  x  —  cst .  With 

ii( y,  t )  =  u(x-  cst,  y,  t )  =  u(x,  y,  t ) ,  the  acceleration  tenn  becomes 

d2u  2  d2u  ~  d2u 

— t  =  G  — t —  2c - 1 — 

dt2  s  d£2  s  d@t  dr 

Similarly  with  v(£,y,t)  =  v(x,  y,  t) 

d2v  2  d2v  ~  d2v  S2v 

- T  —  7 - 2c - 1 - 7T  . 

dt2  s  df-  d@t  dt 2 


square  distortion  speed  \  =  ~ 
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Thus  a  dimensionless  slipper  speed  is  necessary.  To  introduce  this  dimensionless  speed, 


Q 

the  slipper  speed  is  scaled  by  the  distortion  speed  in  region  I  to  define  <x>  =  n2=  —  . 

C\ 

The  last  dimensioned  quantity  which  is  scaled  is  shear  modulus  because  it 
naturally  occurs  in  the  interface  conditions,  2.5.3  and  2.5.4.  Both  the  shear  modulus  in 
region  II  and  region  I  are  scaled  by  the  shear  modulus  in  region  I,  which  produces  the 

Q 

dimensionless  parameter  G  =  =  —  .  Now  the  parameters  values  from  Table  2  are 

Gi 

used  to  derive  the  acceptable  ranges  of  the  dimensionless  parameters.  These  ranges  are 
given  in  Table  3. 


Symbol 

Range 

■S: 

II 

1.5x1  O'3- 1.8x1  O'2 

c  =  cJc2 

0.6-0.94 

a=G,/a, 

33-62 

<o  =  cJc  1 

2-3 

0.3-0.45 

v2 

0.29 

Table  3.  Dimensionless  Parameters  and  Ranges  of  Interest 


The  independent  and  dependent  variables  are  also  scaled.  First  the  £  and  y 
coordinates  are  scaled  by  the  lengths  scales  that  correspond  to  their  direction  of  motion. 

h 

Additionally,  time  is  scaled  by  — ,  the  amount  of  time  it  takes  the  distortion  wave  in 

ci 

region  I  to  travel  through  the  thickness  of  the  coating.  This  implies  the  dimensionless 
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£  y  c 

independent  variables  are  ^  =  — ,  rj  =  — ,  and  z  =  —t .  With  this  transformation,  //  =  0  at 

l  h  h 

the  interface  and  77  =  1  at  the  top  of  the  coating. 

Finally,  the  dependent  variables,  u  and  v ,  are  scaled  in  the  same  way  as  £  and 
.  That  is,  define 


u(x,y,t )  =  u 


f  x  -  cst  y  cxt  ^ 

/  VT 


=  u(£,tj,z) 


2.7.1 


and 


\v{x,y,t)  =  v 
h 


( x-cjt  y  cxt  '' 
I  ’~h’~h 


=  v(C,rj,T) 


2.7.2 


Now  the  moving  coordinates  and  the  set  of  dimensionless  parameters  are  used  to 
transfonn  the  equations  of  motion  and  the  auxiliary  conditions  into  an  equivalent 
nondimensionalized  system.  By  applying  2.7.1  and  2.7.2  to  2.2.1,  2.2.2,  2.2.3,  and  2.2.4, 
the  dimensionless  equations  of  motion  in  the  coating  become 


<5!(2sW9i~<r)  |p  +  g(7)  |p + S'g(il)(n  +  ?, ) 


di<  =  -2&,1“l  +  £A  2.7.3 


d(di  dr 


S2  (g(7)  -  )  |p  +  2gW?,  |p +  SOX'!  +  <?, )  pp 


=  -2c5h> 


fi2v,  ,  d2vt 


-  +  - 


fir" 


2.7.4 


and  in  the  rail  become 


52n,  S2i7 


^•2  (2^2 -c2)— f +  — -f +  ^2(/y  +^r2) 


fi2w 


dC,  drj 


=  -2  —  8 


C  „  5  H,  C  fi“77 


-  +  - 


54^77  co  dt^dz  co2  dz2 


2.7.5 


_2/.  2\dv,  ,  5  v,  .  .  cFw.,  -C  <3  v,  c  5  v, 

^(l-c2)^  +  2<72— K(r2+<72)  — ^  =  -2-£  — ^  +  - 


dij “ 


d^drj  co  d^dz  co1  dz 2 


2.7.6 


where 
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r.  =  - 


1  -  2v, 


;  q,  = 


1-V| 
1  -  2v. 


z  =1,2 


2.7.7 


Now  the  same  transformations  are  applied  to  the  continuity,  boundary,  and  initial 
conditions.  For  continuity  of  displacement,  these  equations  look  exactly  the  same  as  the 
original  functions: 

ux(£,0,t)  =  u2(C,0,t)  2.7.8 


v1(^',0,r)  =  v2(4',0,  r) 


2.7.9 


It  is  now  important  to  write  stress  in  terms  of  2.3.4  to  include  the  scaled  variables  and 
parameters.  Thus  at  77  =  0 ,  the  continuity  of  shear  stress  becomes 


du 


— 1  +  8  —  =  G 


dvx 


du. 


dr/  dC,  \dr)  d£ ) 


and  the  continuity  of  direct  stress  at  77  -  0  becomes 


dux  dvx 

r,  — L  +  qx  — 1  =  G 
dC,  drj 


(  du1  dv,  ^ 

dQ  d  77 


At  the  upper  boundary  of  the  coating,  77  =  1 ,  the  loading  conditions  become 


dA+S1^7=Sh(,r)H{  1-Id) 


and 


di]  dC, 


r'Jc+q'Jr]=  r)// (1  ~  1^ 


2.7.10 


2.7.11 


2.7.12 


2.7.13 
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from  (2.6.1)  and  (2.6.2)  where  / (C,t)  =  — / 

Gi  v 


1  J x-ct  c,t ^ 


s  1  J  and 


I  h 


k(C,r)  = 


1  ,( x-cst  cxt ^ 

/  ’T 


2G, 


i  v 


.  Notice  that  r.  and  q,  are  dimensionless  and  the  functions 


/  and  A:  have  the  same  dimensions  as  Gx ,  thus  /  and  £  are  also  dimensionless  as  their 
strength  has  been  scaled  by  the  shear  modulus  (Gj)  in  the  coating.  The  vanishing  stress 
conditions,  a  . ,  <r .  — >■  0  as  l^l  — >  oo  and  a  — >  0  as  77  — >  -00 ,  are  unchanged  from  the 


transfonn.  Now  at  time  r  =  0  ,  the  scaled  quiescent  initial  conditions  become  u  =  v  =  0 
and 

du  0  du  dv  _ 

- —  = - cm> —  =  0 

dC,  dr  dC, 

2.8  Steady  State  System 

The  steady  state  solution  comes  from  the  system  which  is  void  of  all  time 

dependence.  This  implies  as  time  progresses  the  system  must  see  no  change  in  the 

displacement  with  respect  to  time  or  position  in  the  direction  of  motion.  Mathematically, 

the  right  hand  side  of  the  equations  of  motion  must  be  zero  which  implies 

du ,  ou  ■ 

-2  coS^  +  -  J 


and 


dC,  dr 


dv i  dv, 

-2a>8 — -4 - - 

d£  dr 


must  both  be  constant.  This  thesis  assumes  the  source  of  the  impact  load  moves  at  a 
constant  velocity.  Physically  this  means  the  slipper  has  no  acceleration  because  it  is 
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traveling  at  its  constant  speed.  Without  acceleration,  this  research  assumes  the  system 
has  reached  steady  state  and  the  equations  of  motion  in  the  coating  become 


8 2  ( 2 g(q)ql  -  or )  +  g(rf)  +  82g(q)(q  +  qx )  =  0 


S2  ( g(r/)  ~  or )  +  2g(r/)ql  +  g(r/)(>\  +  qx )  ^  =  0 

v  ’  dC  dr?  d@q 


and  in  the  rail  become 


e2 2 \d2u2  d2u0  ,,2/  x  32v, 


2.8.1 

2.8.2 


2.8.3 

2.8.4 


2.9  Summary 

This  chapter  outlined  the  general  theory  governing  the  problem  of  interest.  The 
theory  led  to  a  system  of  partial  differential  equations  in  terms  of  displacement  to  model 
the  rail  (region  I)  and  the  coating  (region  II)  along  with  the  load.  Then  the  conditions 
necessary  to  solve  the  system  2.5. 1-4  and  2.6. 1,2  were  presented.  Next,  the  parameters  of 
the  system  were  used  to  find  dimensionless  quantities  and  the  parameters  were  scaled  to 
create  dimensionless  parameters  (Table  3).  Finally,  the  steady  state  equations  were 
presented  as  the  governing  equations  2.8. 1-4  to  be  solved  in  this  study. 


24 


Chapter  3.  Methodology  and  Solutions 

3.1  Introduction 

This  chapter  describes  the  methods  used  to  solve  the  system  of  equations 
presented  in  Chapter  2.  For  comparison,  the  problem  is  first  solved  with  the  region  II 
only,  so  that  there  is  no  coating  and  the  force  is  applied  directly  to  the  rail.  Since  the 
system  has  an  infinite  domain  in  the  direction  of  motion  and  vanishing  stress  applies,  a 
Fourier  transfonn  is  used  to  convert  the  equations  to  a  system  of  ordinary  differential 
equations.  Since  the  loading  is  compactly  supported,  the  system  can  be  readily  solved  in 
the  Fourier  domain.  Then  the  solution  technique  is  described  and  the  analytic  solution  is 
presented.  Then  an  inverse  Fourier  transform  is  applied  to  obtain  the  solution  in  the 
spatial  domain. 

The  transformed  equations  for  the  problem  with  coating  are  presented  without 
detail  since  the  method  is  identical  to  the  no  coating  problem.  A  finite  difference  method 
is  introduced  to  solve  the  equations  numerically.  This  method  produces  solutions  to  the 
problem  where  the  shear  modulus  varies  through  the  thickness  of  the  coating.  The  results 
for  variable  coating  thickness  and  shear  modulus  are  compared  via  their  influence  on  the 
stress  distribution  at  the  interface. 

3.2  Rail  Only 

In  this  section,  the  problem  is  solved  with  no  coating  applying  the  boundary 
conditions  directly  to  the  rail.  The  new  boundary  conditions  are  presented.  A  Fourier 
transform  is  introduced  in  the  C,  direction  and  a  homogenous  linear  system  of  first  order 
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differential  equations  is  set  up  and  solved.  Finally,  the  transfonn  is  inverted  and  the 
solutions  for  stress  are  presented  rather  than  displacements. 

3.2.1  New  Boundary  Condition 

The  equations  of  motion  for  the  uncoated  region  II  remain  the  same  and  are  given 
by  2.8.3  and  2.8.4.  However,  with  the  coating  removed,  there  is  no  longer  a  need  for 
continuity  conditions  at  the  interface  of  the  regions.  Instead,  the  force  is  applied  directly 
to  the  half-space.  That  is,  the  upper  boundary  is  at  rj  —  0  ,  and  the  continuity  of  stress 
conditions  at  this  boundary  are 


G 


du 


-CL  +  S 
dr]  df 


Sv,  Sf(() 


ff(Hd) 


3.2.1 


and 


du0  dv, 

— -  +  q2  — 

dC  dr ] 


KQ 


ff(Hd) 


3.2.2 


3.2.2  Fourier  Transform 

To  create  a  system  of  ordinary  differential  equations,  first  a  spatial  Fourier 
transform  is  introduced  in  the  f  direction.  This  transfonn  is  applied  to  the  equations  of 
motion  as  well  as  the  boundary  conditions.  Since  the  rail  is  modeled  as  infinite  in  both 
directions,  the  Fourier  transforms  for  the  functions  u  and  v  are  defined  as 


U(a,j])=  j  eia(u(C,7])dC  3.2.3 


and 


V(a,r/)  =  \  eia(v(C,?])dC  3.2.4 
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Notice,  with  this  transform,  all  variables  and  parameters  except  for  C,  are  considered 
constant.  Further,  since  the  Fourier  transform  is  a  linear  operator,  each  tenn  in  the 
equations  can  be  transformed  separately. 

For  demonstration  purposes,  the  transform  is  perfonned  on  each  tenn  type  in 
2.8.3  and  2.8.4  with  respect  to  the  function  u  : 

]e“ld!“g(p'/W  =  -a2U(a,,!) 

f  jat  ,ij)  j _  d2U{a,rj) 

e  ^  2  ^  ^  2 
_oo  O7)  O7) 

leta(8^,^dc=iaau(a,n) 

^  d^drj  drj 

The  boundary  and  continuity  conditions  are  similarly  transformed. 

Thus,  for  every  derivative  of  u  and  v  with  respect  to  two  derivatives  of  £  ,  the 

transfonned  function  multiplies  -a2  and  its  coefficient.  And  for  every  derivative  of  u 
and  v  with  respect  to  one  derivative  of  ^  ,  the  transformed  function  multiplies  ia  and  its 
coefficient.  This  produces  the  new  system  of  ordinary  differential  equations  in  the 
Fourier  domain: 


-a 


l82  ( 2q2  -c2)U2  +  d-^  +  iaS2(r2  +  q2)d^  =  0 
y  ’  at]  dr] 


-a 


’S2(l-c2)V2  +  2q2  +  ia(r2  +  q2 )  <^=-  =  0 


with  the  following  conditions  to  satisfy: 


•  TT  dV2  A 

iar2U2  +  q2  — -  — >  0 ,  as  77  — >  -00 
dr/ 


3.2.5 

3.2.6 

3.2.7 
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and  at  77  =  0 


+  iaS2V2  =  —F(a) 
dp  G 


■  TT  dV2 

iar2U2  +  q2  — - 
dp 


- —K(a ) 
G 


3.2.8 

3.2.9 


where 


F(a)  =  |  eia(f(C)dC  3.2.10 


and 


1 

K(a)  =  \eiaik(C)dC  3.2.11 

-1 


3.2.3  Eigenvalue-Eigenvector  Approach 

To  find  solutions  for  U  and  V ,  the  equations  of  motion,  3.2.5  and  3.2.6,  are 
converted  to  a  homogeneous  linear  system  of  first  order  differential  equations  in  the 
transfonn  domain.  First,  let 


W  = 


u 

V 

Un 

v„ 


3.2.12 


Flere  the  notation  U  indicates  differentiation  of  U  with  respect  to  77  .Then  from  3.2.5 
and  3.2.6, 
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dW 

drj 


~Un~ 

K 

Vnn 

W3 

-a282  [lq2  -  c2)Wj  -iaS2(r2  +q2)W4 

Z^(  1  ~c-)W2 --i-to(rs  +  q2  )W, 
2q2  lq2 


where  Wk  is  component  k  of  the  vector  W  .  Now  let 


Y  = 


a 


0  0 

0  0 

2S2(2q2-c2)  0 

crJ2 


1 

0 

0 


0 

1 

-iaS2(r2  +q2) 


0 


2q2 


(l-c2)  ~  —  ia(r2+q2) 


2q 


3.2.13 


then  the  equation  becomes 


dW 

dq 


=  YW 


3.2.14 


This  system  has  the  solution 


w=y.b 

7=1 


jXje 


where  A.  is  the  j  eigenvalue  of  Y ,  x  j  is  an  eigenvector  associated  with  the  j 


•  th 


eigenvalue,  and  the  B  ’s  are  unknown  coefficients.  Further,  the  notation  will  mean 


the  k,h  component  of  the  j"'  eigenvector.  Then,  by  the  definition  of  W 


U2(a,Ji)  =  YJBJXj)e 

7=1 


and 


V2(a,q)  =  YJB]x{2)e 


(2L-V? 


7=1 
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Now,  let  px  -  S^Jl-c2  and  p2=S  1  -  — — ,  then  the  eigenvalues  of  7  are  /l,  2  =  +a px 

V  2ci2 


and  Aj  4  =+ap2  with  corresponding  eigenvectors 


•*i,2 


1 

1 

i 

±  — 

S2 

Pi 

and  jc3>4  = 

^3,4 

2 

.  2 

_ ;/7 

JPi 

82  _ 

Now  the  solutions  are  written  as 


U2(a,rj)  =  Bxea™  +  B2eap'n  +  B \ea™  +  B4eapp] 


3.2.15 


and 


V2{a,p)  =  Bx  —  e~a p'n  -B2—  eap'n  +  B3  ^ ea™  -B^±- ea ™ 


S2 


82 


3.2.16 


Pi  Pi 

With  the  form  of  the  solution  known,  it  is  necessary  to  apply  the  boundary 
conditions  to  find  the  unknown  coefficients.  First,  consider  the  vanishing  stress 
condition  as  ij  — »  -oo .  Applying  3.2.15  and  3.2.16  to  3.2.7  implies 


iar2U2  +  q2  ^  =  -ia  [l B \eam  +  B2eappl  ]  +  zi 


la 


_  _Eki 

2  ^2 


V 


82 


J 


~  B3eap"J1  +  B4eap2n~ 


which  must  vanish  as  77  — >  -oo .  This  gives  rise  to  two  cases  to  consider,  a  >  0  and 
a  <  0 .  First  observe,  px  and  p2  are  always  positive  for  the  values  of  the  parameters  this 
study  discusses  and  77  is  always  negative.  Thus,  for  a  >  0  ,  require  Bt  =  B3  =  0  or  the 
stress  would  grow  exponentially  as  77  — »  -00 .  Likewise,  for  a  <  0 ,  require  B2=B4=  0  . 
This  reduces  the  solutions  for  displacement  to 
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3.2.17 


Ufa,rf)  =  /3^m  +  (3^™ 


and 


V2{oc,rj)  = 


Vi 


Win 


Pi 


8l 


sgn(-a) 


3.2.18 


where  sgn(-a)  is  the  sign  of  -a  ,  and  /?,  and  (32  are  unknown  coefficients. 

3.2.4  Determining  the  Coefficients 

To  find  /?,  and  f3: ,  3.2.8  and  3.2.9  are  applied  to  the  general  solution,  3.2.17  and 


3.2.18.  Let  U,  = 


U 

V 


in  region  k  .  Then  for  convenience  the  boundary  conditions  are 


written  in  matrix  form: 


"1  o' 

dU , 

"  0 

iaS2 

— —  + 

_0  q2\ 

dij 

iar2 

0 

8F 

G 

-K 

~G 


From  3.2.17  and  3.2.18, 


where  f3  = 


A 

A 


and 


U2=Z 


J«|  P\V 


0 

a\p2rj 


P 


x  = 


1 


1 


— sgn(-a)  ^fsgn(-a) 
Pi  $ 


3.2.19 


3.2.20 


3.2.21 


Thus  U2(a,0)  =  xP  and 
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dU2 

drj 


(a.O  )  =  x 


a\p\ 

0 


0 

\a\p2 


P 


Substituting  this  into  3.2.19  results  in  the  equation 


G 


1 

0 


0 

q2 


Z 


a\pl 

0 


0 

a\Pi 


iar2 


iaS2 

0 


which  can  be  reduced  to 


3.2.22 


QP  = 


5F_ 

\cc\G 

-iK 

aG 


with  Q  defined  in  Appendix  A.  Solving  for  /?  and  substituting  into  3.2.20  produces 


U2(a,r/)  =  ^z 


Pin 


0  eWP2n 


Q 


A  o 


0  -- 


F 

K 


3.2.23 


which  implies 


U2(a,jj)  = 


1 


Gdet(  Q) 


i( 

Q2$M +($**) 

-(p2eap'n  -QpaM  ) 

/ 

a v  ' 

a 


SQ, 


|2_g|«| Pit  |  P 2  c\a\P2l 


V  p  1 


s 


\a\ 


Pi 


PiQx  e\a\ P2ri 

52 


F 

K 


3.2.24 


2  .  02  A 

V  o  p 

where  Qx  =  — -  and  Q2  =  r, - \ .  Writing  the  solution  in  this  form  breaks  £/,  into 

P\  <>' 


the  coefficients  of  F  and  K  .  Now  the  influence  of  direct  and  shear  stress  on 


displacement  can  be  seen  immediately. 


3.2.5  Inverting  the  Transform 
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The  inverse  Fourier  transform  of  3.2.24  is  applied  only  to  the  coefficients  of  F 
and  K ,  since  the  Fourier  inverted  coefficients  can  be  convolved  with  the  functions  / 

and  k  outside  the  transform  domain  to  obtain  the  solution.  Further,  this  research  is 
interested  in  stress.  Consequently,  the  equations  for  stress  in  the  transform  domain  are 
written  using  the  displacement  solutions  and  the  inverse  transform  is  applied  to  the  stress. 
Similar  to  the  left  hand  sides  of  3.2.8  and  3.2.9,  the  Fourier  transfonned  equations  for 
stress  are 


dCr  =2  G 


iaq2U2  +  r2 


dV2 

drj 


a  =2  G 

nn 


Fn=G 


f  ciV  ^ 

iar2U2  +  q2  — - 
d? 7 


dU2 

dq 


+  iaS~V ' 


J 


Now  substituting  the  solution  3.2.24  into  these  equations  gives 


°K  ~ 


det(Q) 


Q2emm  + 


Pi  " 
\  5  J 


J  a\Pil 


f 


-p2<r™+Qt 


.2  \ 


P2 

q2-^r2 

v  8  J 


J MPin 


iS  sgn(a)  0 

"F" 

0  1 

K 

3.2.25 


o  = 

nn 


det(Q) 


-Qnd 


»|  Pin 


P.2 

8 


-2  \ 


r, - ~-q2 


P2e"m+Ql 


Pi 

82 


J 

2  A 


r> - ~~q2 


J  a\p2n 


,\a\P2n 


iS sgn(a)  0 

~F~ 

0  1 

K 

3.2.26 
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q2P*  +  S~  +  2 pieW™ 

T 

„  1 

Pi 

" 8  o 

~F~ 

T-]  ~  det(Q ) 

Px+S2  e\«\P«  2Q^Pi’1 

0  ip 2  sgn(a) 

K 

Pi 

Notice,  other  than  constants,  there  are  two  term  types  that  occur  in  the  stresses; 
ik]  sgn(a)ea kin  and  kx e1"^2'1 .  Therefore,  the  inverse  Fourier  transform  is  performed  on 

each  term  type  and  the  results  are  applied  to  3.2.25,  3.2.26,  3.2.27. 

Define  the  inverse  Fourier  transform  as 


oo 

=  —  \f(a)e**da 


Then 


oo 

hi 


2  n 

ik] 
2  n 


&  1  (ikx  sgn(a)e“k2,/ )  = 
ikx  sgn(a)ea kl’} e~‘ai da  = 

ea(k2n-i()  +ea(-k2n-i()da  = 


k  < r 

ft  kff+C2 


and 


A. 

2n 


I 


J-\kle^k*)  = 

oo 

—  f  k/^e^da  = 

2ft  J 

-oo 

ea(k2n-i()  +  ea(-k2n-i()da  = 


Kk2  7 
ft  klrf+C2 
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Let  Z,(<r,i?)=  ,  f  , 
pfr+C 


P 


for  i  =  1, 2 ,  then  inverse  Fourier 


and  Nt{^,rf)  = 


pW+C2 


transfonns  for  the  stresses  are 

2 

acc  ~ 


ndet(Q) 


8 


Q2Z1(C,tJ)  + 


r  Pi  a 

<h~^Tr2 
v  s  J 


z2(C,p) 


V+ 


Pi 


pMC,ti)-qx 


Pi 

g, - f  r7 

§2  2 


n2(C,p) 


a nn 


ndet(Q) 


f 

(  2  \ 

\ 

8 

-02zi(C>g)+ 

Pi 

I'l-^Tdl 
l  S  J 

z2(C,p) 

V 

J 

- f  + 


Pi 


-pMC,p)-Qi 


Pi 

r2  -^2<h 

V  5  J 


n2(C,i?) 


Tc,= 


S(Q,(p;+S2)Nl((,^  +  2plNz((,^Yf- 


ndet(Q)  L 
PiQx  (Zi  (C,  P)  ~  2^2  (C,  V))  *  k 


where  (*)  is  the  convolution  operator  defined  as 

oo 

/*«=J7(<r)g(<r-*)Ar 


3.2.28 


3.2.29 


3.2.30 


Recall  that  the  functions  /  and  k  are  the  loading  functions  at  the  surface.  In  this  study, 
three  loading  functions  are  used  to  describe  the  stress  distribution  in  the  rail  in  Chapter  4. 
The  convolutions  of  the  solution  in  the  rail  with  the  loading  functions  are  given  in 
Appendix  B. 
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3.3  Coated  Rail  Analytic  Solution 

In  this  section,  the  problem  is  solved  with  the  coating  on  the  half-space.  The 
solution  procedure  follows  that  of  section  3.2.  The  solution  in  the  spatial  Fourier  domain 
is  introduced  and  a  homogenous  linear  system  of  first  order  differential  equations  is  set 
up  and  solved  for  the  rail  and  coating.  Unknown  coefficients  are  found  by  applying  the 
boundary  and  interface  conditions. 

3.3.1  Solution  in  the  Fourier  Domain 

To  write  the  equations  of  motion  in  region  I  in  matrix  form,  first  divide  2.8.1  and 

2.8.2  through  by  g(r/)  and  make  some  convenient  substitutions.  Let 

(  2  x 

a  =  a(rj)  =  a282  2 q - 

l  gin)) 

bx  =  iaS2(r  +  q) 
b2  =  ia(r  +  q ) 

(  2 

c  =  c(rj)  =  o28 2  1  — — — 

l  g(v)j 

where  co  and  8  are  defined  in  Table  3  and  r  and  q  are  define  in  2.7.7.  As  with  the 
uncoated  problem  equations,  2.8.1  and  2.8.2,  are  spatially  Fourier  transformed  to  produce 

for  0  <  77  <  1 .  The  boundary  conditions  at  the  top  of  the  coating,  77  =  1 ,  is  similar  to 
3.2.19  and  becomes 


1  0 1  d2U,  To  bAdU, 

0  2  q  dr)2  b2  0  drj 
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where  F  and  K  are  defined  in  3.2.10  and  3.2.1 1.  At  //  =  0 ,  the  continuity  of 
displacement  becomes 

U i(«,0)  =  U2(a,  0) 
and  the  continuity  of  stress  becomes 

T  n  i  rvfid  1 

U1(a,0)  = 


3.3.3 


"1  o' 

1  (a,0)  + 

"  0 

iaS2 

0  q,_ 

dr/ 

iar | 

0 

/ 

"1  o' 

dU , 

0 

iaS 2 

A 

G 

V 

0  t/2_ 

— -(«,0)  + 
dr/ 

iar2 

0 

t/,(«,0) 

7 

3.3.4 


This  system  of  equations  is  converted  to  a  linear  system  of  differential  equations  using 
3.2.12  to  produce  3.2.14  with  7  now  defined  as 


7  = 


f 


2  c-2 

a  o 


co “ 


,2  A 


0 
0 

2 q.  - 

v  s(v)j 

0 


a28 2 


0 

0 

0 

1- 


co~ 


m 

2*i 


1 

0 

0 


0 

1 

iaiy(i]  +qx) 

0 


3.3.5 


2*,  V  i(7)y 

If  there  is  a  variation  in  the  shear  modulus  ( g(r])  is  non  constant),  then  the  system  must 
be  solved  numerically.  However,  if  the  coating  is  uniform  then  g(//)  =  1 ,  and  the  system 
has  a  general  solution  similar  to  the  no  coating  problem.  Thus  the  general  solutions  for 
displacements  in  the  coating  are 


U1(a,r/)  =  YJAJXje 

2=1 


3.3.6 


and 


V\  {a,Tj)  =  YJAjxye 


I2 


3.3.7 


2=1 
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Let  mx  =  iS\lof  -1  and  m2  =  8  1-——  ,  then  the  eigenvalues  are  \2=  +aml  and 

V  2(h 

A3  4  =  +am2  with  associated  eigenvectors 


1 

1 

1 

+  im2 

± — 

and  x, .  = 

8 2 

m. 

1 

K.2 

3,4 

^3,4 

.  2 

im2 

—ia 

■ - 

L  - 

82  _ 

The  general  solution  for  the  rail  is  given  by  equations  3.2.17  and  3.2.18.  All  that  remains 
is  to  find  A1,A2,A2,A4,j3l,  and  J32 . 

3.3.2  Determining  the  Coefficients 

The  coefficients  Al,A2,A3,A4,(3l,  and  /?,  are  determined  by  satisfying  the  load 

conditions  at  the  top  of  the  coating,  3.3.2,  and  the  interface  conditions,  3.3.3  and  3.3.4. 
After  these  equations  are  transformed  and  the  general  solutions,  given  by  3.2.17,  3.2.18, 
3.3.6,  and  3.3.7,  are  introduced,  a  linear  system  for  the  coefficients  emerge.  This 
equation  has  the  form 


R2  ' 

A 

~0~ 

L 

®2x2_ 

F 

Here  A  is  a  vector  of  the  four  unknown  coefficients  associated  with  the  coated  region 
solution  and  j3  the  coefficients  associated  with  the  rail  region  solution.  The  first  four 

\8F~ 

elements  of  the  forcing  vector  are  all  zero  while  the  last  are  comprised  in  F  =  , 

-K 

where  F  and  K  are  defined  in  3.2.10  and  3.2.1 1.  The  matrix  is  comprised  of  sub 
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blocks,  Rx,  R2,  L  given  in  Appendix  C.  Rx  and  R2  together  fonn  a  4x6  block  derived 


from  the  interface  conditions  at  77  =  0 .  The  L  block  is  derived  from  the  loading  stress 
conditions  at  77  =  1 . 

With  the  aid  of  Mathematica™,  the  solutions  to  3.3.9  can  be  determined. 
However,  the  results  are  cumbersome  and  not  suitable  for  transform  inversion.  Instead 
the  solution  to  3.3.9  is  found  for  selected  parameters.  The  resulting  coefficients  can  be 
used  in  3.2.17,  3.2.18,  3.3.6,  and  3.3.7  to  determine  the  displacements  (in  the  transform 
domain)  through  the  thickness  of  the  coating  and  in  the  rail.  These  results  are  used  to 
validate  a  numerical  approach. 

3.4  Finite  Difference  Method/Coated  System 

To  handle  coatings  which  may  vary  through  the  thickness,  it  is  necessary  to  use  a 
numerical  method  to  investigate  displacements  in  the  coating  and  their  influence  on  the 
stress  distribution  in  the  half-space. 

3. 4. 1  Method  Derivation 

The  differential  equations  in  the  coating  are  written  in  terms  of  central  differences 
to  obtain  0{h2\  approximations,  where  the  interval  77  =  [0, l]  is  partitioned  into  n 


subintervals  of  size  h  =  1/ .  For  the  rest  of  this  section,  subscripts  are  no  longer  used  to 

identify  region.  Unless  otherwise  noted,  all  constants  and  functions  are  for  the  coating 
(region  I).  First,  let  Uj  =  U  ( 77, )  be  the  displacement  vector  at  the  ith  partition  starting  at 


77  =  0  and  rji  =kh  ,  k  -  0,1,...,  77 ,  be  the  size  of  each  partition.  Then  central  differences 
for  each  derivative  term  are 
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d2U 


drf 


n=*h 


-^{UttI-2U,  +  17,., 


and 


dU 

drj 


ri=rii 


Tdv"-u'-‘ 


Then  from  3.3.1,  the  central  difference  formula  for  the  equations  of  motion  become 
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3.4.1 


for  k  =  0,1,2,..., n  where  k  =  0  corresponds  to  the  interface  between  the  coating  and  the 
rail,  77  =  0  ,  and  k  =  n  corresponds  to  the  top  of  the  coating,  77  =  1 . 

When  z'  =  0,  equation  3.4.1  requires  knowledge  of  U  x  (i.e.  U (-/?) ).  Also  at 
k  =  77 , 3 .4. 1  requires  Un+l  (U  (/  +  /?  )) .  These  values  are  from  ghost  points  used  in  the 

calculations  and  are  determined  from  the  boundary  and  interface  conditions. 

Now  the  boundary  conditions  are  converted  to  central  differences  as  well.  First, 
at  77  =  1 ,  the  condition  becomes 
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Substituting  into  3.4.1  for  to  k  =  n  produces 
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3.4.3 


At  77  =  0 ,  the  boundary  condition  is  more  complex  to  derive  since  the  solution  in 
the  rail  for  the  coupled  system  is  not  known.  However,  the  form  of  the  solution  is  the 
general  form  of  the  solution  in  the  rail  is  given  by  3.2.20.  To  evaluate  the  solution 
numerically  in  the  coating,  it  is  necessary  to  eliminate  dependencies  on  U2(a, 0)  in  the 
rail  using  the  known  form  of  the  rail  solution.  First,  applying  the  boundary  condition 
3.3.3  and  the  known  solution  3.2.20  produces 

Uj(a,0)  =  xP  =  U2(a,0) 


with  x  given  by  3.2.21.  For  the  parameter  ranges  of  interest,  x  is  always  invertible, 
which  implies 

P  =  X-lU1{a,  0)  3.4.4 

Next,  apply  the  continuity  of  stress  condition  in  3.3.4.  First  observe  that  the  right 
hand  side  of  3.3.4  is  simply  the  left  hand  side  of  3.2.22.  Using  3.4.4  for  ft  in  the  left 
hand  side  of  3.2.22  and  substituting  into  3.3.4  produces  the  new  form  of  the  boundary 
condition: 

Ol 

3.4.5 
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dU , 

0  qx 

2  (a,  0) +  //[/,  (a,  0)  = 
di) 
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where  H  captures  the  influence  of  the  rail  on  the  coating.  H  is  given  in  Appendix  D. 


Now,  introduce  the  finite  differences  at  77  =  0  to  produce 

0 
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Substituting  this  into  3.4.1,  with  i  =  0  produces  the  effective  boundary  condition 
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with  G  given  in  Appendix  D. 


3.4.1,  3.4.3,  and  3.4.6  are  the  equations  necessary  for  the  finite  difference  code. 

3.4.1  Code  Validation 

When  6(77)  =  1 ,  the  finite  difference  code  solves  the  same  system  of  equations  as 
the  analytic  solution  described  in  section  3.3,  the  coated  rail  problem.  Thus,  the  finite 
difference  code  can  be  compared  to  the  analytic  solution  for  selected  values  of 
parameters.  The  solutions  are  compared  over  the  entire  range  of  parameter  values  of 
interest  to  this  study.  Figure  4  illustrates  the  convergence  of  the  finite  difference 
algorithm  as  the  number  of  partitions  in  the  coating  increase,  with  parameters  S  ,  ,  v2 , 

c,  co,  G  ,  F  ,  and  K  set  to  1CT3 ,  0.4  ,  0.29  ,  0.9  ,  3 ,  60  ,  1,  and  1,  respectively.  At 
n  =  100  ,  the  curve  associated  with  the  exact  solution  and  that  of  the  finite  difference 
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solution  agree  almost  exactly  over  the  range  of-10<a<10.  These  curves  are  a  plot  of 
the  transformed  U  displacement  in  the  complex  U  plane.  Curves  corresponding  to 
n-  5,1 0,  and  1 00  are  shown  for  comparison. 
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Figure  4.  Demonstration  of  Finite  Difference  Code 
Convergence  for  Select  Parameter  Values 

During  the  test  for  convergence,  it  was  discovered  that  a  certain  parameter 
combination,  the  product  a 5 ,  had  the  greatest  affect  on  error.  As  the  product  ad  grows 
in  magnitude,  the  number  of  partitions  of  the  interval  0  <  77  <  1  needed  to  maintain  a 
specified  accuracy  must  grow  as  well.  This  study  seeks  to  have  the  relative  error  less 
than  10  4 .  Figure  5  shows  the  number  of  partitions  needed  for  convergence  to  the  true 
solution  within  a  specified  level  of  accuracy  versus  the  product  aS .  The  level  of 
accuracy,  given  by  relative  error,  is  labeled  on  each  curve.  For  example,  with  aS  =  0.7  , 
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follow  the  dashed  line  from  the  aS  axis  up  to  the  curve  corresponding  to  a  relative  error 
of  10  5 .  Then  follow  the  dashed  line  to  the  axis  corresponding  to  the  number  of 
partitions  to  get  n  =  26  ,  partitions  needed  for  ad  =  0.7  with  a  relative  error  of  10  5 . 
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Figure  5.  Number  of  Partitions  Needed  for  Code  Convergence 
With  Respect  to  Relative  Error 

For  this  study,  n  =  15  is  an  acceptable  number  of  partitions  to  obtain  relative  error  less 
than  10  4 . 

3.5  Summary 

This  chapter  began  with  a  discussion  of  the  rail  only  problem.  A  spatial  Fourier 
transform  was  introduced  in  the  C,  direction  which  led  to  a  homogeneous  linear  system 
of  ordinary  differential  equations  with  a  known  general  solution.  Next,  the  coefficients 
for  the  general  solution  were  determined  and  the  solution  was  written  in  terms  of  stress 
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giving  equations  3.2.25-27 .  Next  the  Fourier  transform  was  inverted  to  produce  3.2.28- 
30  in  tenns  of  stress. 

After  the  rail  only  solution  was  presented,  the  same  methods  were  used  to  solve 
the  coating  problem,  which  led  to  the  system,  3.3.9,  to  solve.  Next  the  finite  difference 
method  was  introduced  as  a  method  to  calculate  the  stresses  when  the  shear  modulus 
varied  through  the  thickness  of  the  coating.  The  finite  difference  approximation,  3.4.1, 
for  the  equations  of  motion  and  the  boundary  conditions,  3.4.3  and  3.4.6,  were  derived. 
Finally,  the  code  was  validated  by  comparison  with  the  analytic  solution  in  the  transfonn 
domain.  Then  the  number  of  partitions  of  the  coating  thickness,  n  =  15 ,  for  the  finite 
difference  method  to  converge  to  a  minimum  relative  error  of  10  4  was  demonstrated. 
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Chapter  4.  Results  and  Analysis 


4.1  Introduction 

This  chapter  investigates  the  nature  of  stress  developed  in  the  coating  and  the  rail. 
First  the  problem  with  no  coating  is  examined  under  the  influence  of  several  stress 
distributions  in  the  spatial  domain.  A  point  source,  a  uniform  distribution  over  the  length 
of  the  shoe,  and  a  parabolic  distribution  of  stress  are  studied  to  understand  how  these 
distributions  affect  stress  in  the  uncoated  half-space.  This  is  accomplished  by  using  the 
analytic  solution  for  the  rail  only  problem. 

The  parabolic  loading  distribution  is  then  used  on  the  coated  rail.  Since  the 
solution  for  the  coated  problem  is  numeric,  a  Fast  Fourier  Transform  is  used  to  invert  the 
transfonn  and  study  the  results  in  the  spatial  domain  for  different  constant  shear  moduli 
and  various  coating  thickness.  Then,  using  the  finite  difference  code,  the  shear  modulus 
is  allowed  to  vary  as  a  continuous  function  through  the  thickness  of  the  coating.  This 
continuous  function  allows  the  shear  modulus  to  nearly  vanish  at  the  surface  of  the 
coating,  simulating  a  liquid-like  interface.  It  is  also  used  to  investigate  a  layered  coating 
effect. 

4.2  No  Coating  Results 

4.2.1  Delta  Distribution 

For  a  point  source,  the  force  is  applied  to  a  single  point  on  the  surface  of  the  rail, 
C,  -  0 ,  with  the  use  of  the  Dirac  delta  function  for  the  loading.  The  shear  loading  is 
introduced  through  the  frictional  effects.  Thus,  the  shear  load  is  taken  to  be  a  fraction  of 
the  direct  load  based  on  a  coefficient  of  friction  giving  /  =  0.2 S(£)  and  k  =  in 
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equations  2.7.12  and  2.7.13.  Equations  3.2.28,  3.2.29,  and  3.2.30  are  used  to  produce  the 
following  contour  plots  for  the  point  source  in  the  spatial  domain. 


Figure  6.  Point  Source  Direct  Stress  ( crft„)  Contours 


In  Figure  6,  the  contours  represent  direct  stress,  ,  in  the  direction  of  motion  of 
the  slipper.  The  horizontal  axis  is  the  direction  of  motion,  ^  ,  with  the  positive  to  the 
right  while  the  vertical  axis  is  through  the  thickness  of  the  rail,  77 ,  with  the  negative 
down.  Even  though  the  coating  is  not  present,  the  no  coating  problem  is  still  scaled 
relative  to  a  coating  thickness.  The  graph  shows  a  width  of  one  shoe  length  and  a  depth 
of  one  times  the  thickness  of  a  coating  assuming  5  =  hl  l  —  1 .8x10  2 .  The  labels  on  each 
contour  are  the  relative  magnitudes  (recall  stress  is  scaled  by  the  shear  modulus)  of  the 
stress  for  that  curve.  Notice  in  front  of  the  source,  the  relative  magnitudes  of  the  stresses 
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are  greater  per  depth  into  the  rail  than  that  of  the  stresses  behind  the  source.  Further  the 
sign  of  the  stress  in  front  is  negative  while  the  sign  is  positive  behind.  This  signifies  the 
stress  in  front  of  the  slipper  is  in  compression  while  the  stress  behind  the  slipper  is  in 
tension.  As  C,  approaches  zero,  the  magnitude  of  the  stress  becomes  infinite  because  of 
the  singularity  at  the  origin.  This  is  of  course  due  to  the  singular  nature  of  the  loading 
function.  Further  as  C,  moves  away  from  the  origin,  the  stress  rates  diminish  rapidly. 
This  is  because  the  point  source  concentrates  the  applied  force  at  the  origin. 


Figure  7.  Point  Source  Direct  Stress  ( cr ;)/ )  Contours 


Figure  7  illustrates  the  direct  stress,  <7l/n  ,  through  the  thickness  of  the  rail  in  the 

spatial  domain.  The  strength  of  the  stress  per  depth  is  greater  in  front  of  the  source  just 
as  with  the  stress  in  the  direction  of  motion,  cr  „  .  The  sign  of  the  stress  is  negative  in 
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front  of  the  slipper  showing  a  compression  through  the  thickness  of  the  rail.  The  sign  is 
negative  behind  the  slipper  demonstrating  tension.  Just  as  with  the  direct  stress  in  the 
direction  of  motion,  the  spike  in  the  center  signifies  an  increasing  stress  as  ^  approaches 
zero  due  to  the  singular  nature  of  the  loading  function. 


Figure  8.  Point  Source  Shear  Stress  ( )  Contours 


Figure  8  illustrates  the  shear  stress,  t-  ,  in  the  spatial  domain.  Notice  the 

magnitudes  of  the  shear  stress  per  depth  of  penetration  are  significantly  smaller  than  the 
direct  stresses.  This  is  because  the  shear  loading  was  considered  to  be  a  fraction  of  the 
direct  load.  Additionally,  the  sign  of  the  shear  stress  in  front  of  the  slipper  is  negative 
indicating  the  shear  stress  is  in  compression.  At  the  same  time,  the  sign  of  the  stress 
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behind  the  slipper  is  positive  indicating  the  direction  of  stress  is  in  the  direction  of 
motion.  Further,  the  unlabeled  contours  in  the  center  of  Figure  8  have  the  same 
magnitude  with  opposite  signs  as  the  contours  from  the  top  on  the  same  side  of  the 
source.  This  is  the  behavior  of  the  material  in  small  neighborhood  near  f  —  0 ,  which  is 
directly  beneath  the  point  source.  Just  as  with  the  direct  stresses,  the  shear  stress 
experiences  the  same  singularity  at  the  origin. 

4.2.2  Uniform  Distribution 

With  the  uniform  distribution,  a  constant  load  of  magnitude  one  is  applied  to  the 
entire  length  (0.2  meters)  of  the  slipper  so  that  /  =  0.2//(l  -  l^l)  and  k  =  H(  1  -  \f\) . 


Figure  9.  Unifonn  Source  Direct  Stress  (cr~)  Contours 


Figure  9  illustrates  the  behavior  of  the  direct  stress,  cr  „  ,  in  the  direction  of 
motion  where  f  ranges  over  two  times  the  length  of  the  shoe  and  the  depth,  ij ,  is  100 
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times  the  coating  thickness.  Since  the  uniform  source  is  not  as  concentrated  as  the  point 
source,  it  is  necessary  to  look  deeper  into  the  rail  to  understand  the  behavior.  The 
spreading  of  force  from  the  uniform  source  enables  the  stress  to  distribute  over  more  of 
the  rail  than  the  point  source  as  well  as  limiting  the  penetration  through  the  thickness. 
The  value  of  the  stress  is  negative  from  the  back  of  the  slipper,  ^  =  -1 ,  through  the 
positive  C,  axis.  This  implies  the  direct  stress  is  in  compression  under  and  in  front  of  the 
slipper.  The  stress  changes  sign  at  the  rear  of  the  slipper  implying  the  stress  is  in  tension 
behind  the  slipper.  Also,  notice  that  the  magnitude  of  the  stress  is  greatest  propagating 
from  beneath  the  slipper,  -1  <  C,  <  1 .  Figure  10  illustrates  the  same  behavior  for  the 
direct  stress  cr;/;;  as  a  .  Since  the  downward  direction  is  the  negative  77  axis,  this  also 

implies  the  stress  is  in  compression  beneath  and  in  front  of  the  slipper  while  the  stress  is 
in  tension  behind  the  slipper. 


y 
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Figure  10.  Uniform  Source  Direct  Stress  Contours 


Figure  11.  Uniform  Source  Shear  Stress  (zv;J  Contours 


Figure  1 1  illustrates  the  shear  stress,  t~  ,  with  a  unifonn  source.  Notice  that  if 

the  load  were  not  moving,  the  graph  would  be  perfectly  symmetric  and  the  stress  would 
be  zero  at  the  center  of  the  slipper,  C,  =  0 .  Since  the  load  is  moving,  the  location  of  zero 
stress  is  shifted  in  the  positive  direction  of  motion.  In  front  of  the  slipper,  the  sign  of  the 
stress  values  are  negative  indicating  the  shear  stress  is  in  compression.  Behind  the 
slipper,  the  sign  is  positive  indicating  the  shear  stress  is  in  tension.  Further,  behind  the 
slipper,  the  stresses  of  equal  magnitudes  as  in  front  propagate  further  through  the 
thickness,  down  the  negative  77  axis. 

4.2.3  Parabolic  Distribution 
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With  the  parabolic  distribution,  the  load  is  not  constant  over  the  length  of  the 
slipper.  The  load  is  greatest  at  the  center  of  the  slipper  and  weakest  at  the  edges.  Define 
the  formula  for  the  parabolic  distribution  as 


f,.(x) 


—  (1~) 
4  L  Z2 


r  I  <  L 

Ixl  >  L 


4.2.1 


where  L  is  the  length  of  the  slipper.  The  plots  for  the  parabolic  distribution  use  the 
value  L  =  1 . 


Figure  12.  Parabolic  Source  Direct  Stress  (cr^)  Contours 


Figure  12  illustrates  the  direct  stress,  oct- ,  in  the  direction  of  motion.  All  the 
magnitudes  of  stress  are  negative  indicating  the  stress  is  in  compression.  Further,  as  ^ 
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and  ij  approach  zero,  the  curves  are  more  concentrated.  This  is  because  the  force  due  to 


loading  is  greatest  in  the  center  and  the  strength  of  the  stress  is  weakening  through  the 
thickness.  Figure  13  illustrates  the  direct  stress,  cr;?7  ,  through  the  depth  of  the  rail.  The 

stress  is  of  lower  magnitude  because  it  is  primarily  due  to  the  direct  stress  loading 
function  and  does  not  reflect  both  boundary  conditions  like  cr.-, ,  but  has  the  same 

behavior  as  the  stress  in  the  direction  of  motion. 


Figure  13.  Parabolic  Source  Direct  Stress  ( <Jm )  Contours 


54 


Figure  14.  Parabolic  Source  Direct  Stress  j  Contours 


Figure  14  illustrates  the  shear  stress  in  the  rail  for  the  parabolic  distribution.  One 
interesting  aspect  of  this  graph  is  the  similarities  with  the  shear  stress  for  the  uniform 
source  (Figure  11).  Just  like  the  unifonn  source,  the  motion  shifts  the  graph  to  the  right 
for  the  parabolic  source.  The  value  of  the  stress  is  negative  indicating  the  stress  is  in 
compression. 

The  parabolic  distribution  is  an  interest  in  this  study  because  the  properties  of  this 
distribution  can  be  easily  observed  in  the  problem  with  coating.  When  the  coating  is 
applied  to  the  rail,  it  acts  as  a  transfer  function  modifying  the  distribution  of  stress 
applied  to  the  surface  to  a  new  distribution  applied  to  the  rail  at  ij  —  0 .  Before  the 
coating  problem  is  examined,  first  the  properties  of  the  parabolic  distribution  function  are 
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L 

examined.  Notice  from  4.2.1,  that  if  fL (x)  is  stress,  then  J  fL(x)dx,  or  the  area  under 

-L 

the  curve,  is  force.  The  area  under  this  curve  is  always  one  despite  the  value  for  L  ,  that 
is  the  total  applied  force  is  always  one.  The  area  denotes  the  strength  of  the  applied  load 
distributed  over  a  region  from  —L  to  L  . 

By  increasing  the  value  for  L  ,  the  distribution  of  stress  is  applied  to  a  larger 
surface  without  increasing  the  total  applied  force.  This  implies  a  weaker  distribution  of 
stress  on  the  surface  of  the  rail.  Figure  15  illustrates  the  effects  of  increasing  L  .  As  L 
increases  to  10,  the  stress  of  magnitude  1  is  applied  to  the  length  from  -10  to  10  rather 
than  from  - 1  to  1  when  L  =  1 . 
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Figure  15.  Comparison  of  the  Parabolic  Distribution 
Function  with  Different  L  values 

This  spreading  effect  in  the  stress  applied  to  surface  is  the  expectation  for  the  coating 
acting  as  a  transfer  function. 
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4.3  Rail  with  Coating 


This  section  investigates  the  effects  on  stress  of  coating  the  rail  with  a  thin 
polymer  or  epoxy.  First  the  effects  of  various  thicknesses  and  shear  moduli  are  studied. 
Then  the  distribution  of  stress  is  studied  with  the  shear  modulus  varying  as  a  function  of 
thickness  to  investigate  the  effects  of  a  vanishing  shear  modulus  near  the  surface  of  the 
coating  and  layering  coatings.  The  loading  function  used  is  the  parabolic  distribution 
with  L  =  1 . 


4.3.1  Effects  of  Coating  Thickness  and  Shear  Modulus 

Recall  that  the  coating  thickness  is  reflected  in  the  dimensionless  parameter 
8  =  hi l ,  whose  ranges  are  given  in  Table  3.  The  coating  thickness  is  examined  at  three 
levels,  1.5xl0“3(low),  9. 8x1  CT3  (mid),  and  1. 8x1  CT2  (high).  The  influence  of  the  coating 
on  the  stress  distribution  in  the  rail  will  be  depicted  by  the  stress  distribution  at  the 
interface,  ?j  -  0 .  Consequently,  the  graph  of  the  parabolic  distribution  applied  to  the 
surface  in  the  rail  only  problem  is  plotted  with  the  stress  curves  for  cr  w  and  z>  to  show 

the  differences  between  the  uncoated  and  coated  problems.  The  parabolic  distribution  is 
not  plotted  with  cr~  because  the  boundary  conditions  do  not  include  direct  stress  in  the 

direction  of  motion. 
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Figure  16.  Comparison  of  Direct  Stress  (cr~)  with  Respect  to  Coating  Thickness 

In  Figure  16,  the  horizontal,  C,  ,  axis  is  the  direction  of  motion  of  the  slipper  and 
the  vertical  axis  is  the  direct  stress,  o.c  ,  in  the  direction  of  motion  at  ij  —  0 .  All  stress 

generated  comes  from  the  two  loading  functions,  /  and  k  .  Even  though  the  stress  is 
plotted  with  respect  to  all  three  coating  thicknesses;  there  is  no  distinct  difference 
between  them.  Further,  the  load  was  applied  from  -1  to  1  on  the  £  axis.  All  three 
curves  also  intersect  the  horizontal  axis  at  C,  =  —  1  and  C,  —  1  as  does  the  load.  This 
implies  there  is  no  difference  in  the  spreading  of  the  load  between  the  three  coating 
thicknesses. 
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Figure  18.  Comparison  of  Direct  Stress  with  Respect  to  Coating  Thickness 
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Figure  17  illustrates  the  effects  of  the  coating  thickness  on  direct  stress,  cr;/7  , 

through  the  thickness.  The  load  is  graphed  along  with  the  stress.  Figure  18  illustrates  the 
effects  of  the  coating  thickness  on  shear  stress,  z>  .  Notice  in  both  figures,  the  peaks  of 

the  graphs  do  not  show  a  significant  difference,  so  the  intensity  or  strength  of  the  load 
has  not  changed  due  to  the  coating.  Like  the  direct  stress  in  the  direction  of  motion,  there 
is  also  no  spreading  effect  due  to  the  coating. 

Next  a  constant  shear  modulus  is  varied  to  investigate  its  effect  on  the  stress 
propagated  into  the  rail.  Recall  the  dimensionless  parameter  G  =  G2/Gl  whose  ranges 
are  given  in  Table  3.  This  relationship  implies  that  the  dimensionless  parameter  G 
decreases  as  the  shear  modulus  of  the  coating  increases,  so  the  three  chosen  parameters 
are  62  (low),  47  (mid),  and  33  (high).  In  the  following  plots  the  parameter  5  is  always 
chosen  as  1.5x10  3 .  The  curves  for  G  =  33  are  the  same  curves  as  the  8  =  1 .5x1 0  '  in  the 
figures  for  coating  thickness. 
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Figure  19.  Comparison  of  Direct  Stress  (ov  j  with  Respect  to  Shear  Modulus 


As  illustrated  in  Figure  19,  the  direct  stress,  cr  „  ,  in  the  direction  of  motion  is 

significantly  affected  in  magnitude  by  the  shear  modulus  of  the  coating.  The  area  under 
the  curve  for  G  -  33  or  the  load  is  1.09  while  the  area  under  the  curve  for  G  =  62  is 
0.89.  Recall,  the  area  under  the  curve  is  the  total  application  of  force.  This  is  an 
improvement  of  0.2  in  scaled  force  applied  to  the  rail  at  the  interface.  As  the 
dimensionless  shear  modulus  decreases  (increases  in  the  coating),  the  strength  of  the 
stress  at  the  interface  between  the  rail  and  the  coating  increases.  This  implies  that  a 
lower  shear  modulus  for  the  coating  would  allow  less  direct  stress  in  the  direction  of 
motion  at  the  interface  with  the  rail.  Because  the  load  is  not  applied  to  the  rail  through 
the  direct  stress  in  the  direction  of  motion,  <t  _,  ,  this  does  not  imply  that  the  stress 

distribution  in  the  rail  has  been  weakened.  However,  this  reduction  in  stress  can  be 
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related  to  material  yield  or  failure.  Through  theory  on  design  stress  or  Von  Mises  stress, 
it  can  be  shown  that  a  reduction  in  the  direct  stress,  <j  x  ,  can  reduce  the  likelihood  of 

material  failure  (20). 

The  figures  for  direct  stress,  a  ,  through  the  thickness  and  shear  stress,  z>  ,  are 

not  given  separately  because  they  are  exactly  the  same  as  Figures  17  and  18,  respectively. 
This  demonstrates  that  the  value  of  the  shear  modulus  in  the  coating  has  very  little  effect 
on  the  direct  stress  through  the  thickness  or  the  shear  stress.  It  was  expected  that  the 
higher  shear  modulus  in  the  coating  would  disperse  the  load  across  a  larger  area  of  the 
rail  due  to  the  faster  stress  wave  speed.  However,  the  coating  is  too  thin  to  see  this 
effect. 

If  the  shear  modulus  in  the  coating  starts  to  approach  zero,  or  become  more  liquid 
like,  there  is  a  greater  effect  in  the  direct  stress  in  the  direction  of  motion.  Figure  20 
illustrates  this  behavior  with  G  =  1000 ,  or  Gl  =  8x10  2  kg/ms2,  which  is  an  unrealistic 
value  for  a  polymer  or  epoxy  adhesive. 
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Figure  20.  Comparison  of  Direct  Stress  {<ycA  with  an 
Extremely  Low  Shear  Modulus  in  the  Coating 

4.3.2  Variable  Shear  Modulus  through  the  Thickness 

This  section  investigates  the  affect  varying  the  shear  modulus  as  a  function  of 
thickness  has  on  the  stress  distribution  in  the  rail.  The  function  g(?7 )  in  2.8.1  and  2.8.2 
is  used  to  vary  the  shear  modulus.  A  function  similar  is  used  for  functionally  graded 
materials  (21:615).  Define  g(?7)  as 

g(<l)  =  (0,-Gb)ri’  +G„ 

where  Gh  is  the  starting  shear  modulus  at  the  interface  ,  Gt  is  the  ending  shear  modulus  at 
the  surface  of  the  coating,  and  m  dictates  the  profile  of  variation  through  the  thickness. 
Figure  21  illustrates  the  behavior  of  the  function  g(//)  for  different  values  of  m.  When 
m=  1,  the  coating  properties  vary  linearly  through  the  thickness.  When  m  is  greater  than 
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one,  the  material  maintains  a  shear  modulus  close  to  its  ending  shear  modulus  through  the 
thickness,  and  then  the  modulus  vanishes  close  to  the  surface.  When  m  is  less  than  one, 
the  materials  shear  modulus  vanishes  quickly  as  77  approaches  the  surface  of  the  coating. 


Figure  21.  Behavior  of  Shear  Modulus  Varying  Function  g  ( // ) 

For  Different  m  Values 

The  graphs  using  the  variable  shear  modulus  function  plotted  for  m= 0.5  and  m= 2 
are  illustrated  in  Figures  22,  23,  and  24.  The  best  case  for  the  constant  shear  modulus, 

G  =  62 ,  is  also  plotted  for  comparison.  The  input  for  ^(77)  is  Gt  =  0.001  and  Gh=  62. 
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Figure  24.  Comparison  of  Direct  Stress  for  Variable  Shear  Modulus 

In  each  figure,  the  curves  are  nearly  indistinguishable.  This  illustrates  there  is  no 
significant  difference  in  layering  coatings  with  different  material  properties  versus  having 
one  material  with  a  low  shear  modulus. 

4.4  Summary 

First  the  effects  of  the  point  source,  uniform  source,  and  parabolic  source  on 
stress  distribution  in  the  rail  without  coating  were  demonstrated.  Then  the  effects  of 
coatings  on  the  parabolic  stress  distribution  were  studied.  It  was  discovered  that  the 
presence  of  a  thin  coating,  despite  shear  modulus  or  thickness,  does  not  affect  the  spread 
of  the  distribution  of  force  to  the  rail.  Further,  the  thickness  of  the  coating  had  little 
effect  on  the  magnitude  of  the  force  propagating  to  the  rail.  Also,  it  was  discovered  that 
a  lower  shear  modulus  in  the  coating  material  produces  a  significantly  lower  direct  stress 


66 


in  the  direction  of  motion.  This  does  not  imply  a  weaker  load  is  transferred  to  the  rail 
through  the  boundary  conditions.  However,  it  does  imply  that  the  probability  of  failure  is 
reduced.  This  lower  shear  modulus  implies  the  coating  less  rigid  or  stiff.  Finally,  it  was 
shown  that  layering  coatings  produces  no  benefit  greater  than  that  of  a  constant  shear 
modulus  through  the  thickness  of  the  coating. 
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Chapter  5.  Summary  and  Conclusions 

5.1  Summary 

The  first  chapter  of  this  thesis  discussed  the  motivation  for  this  research,  the 
HHSTT’s  interest  in  mitigating  the  gouging  phenomenon.  Additionally,  the  chapter 
reviewed  past  research  in  half-space  stress  propagation  and  problem  fonnulation,  much 
of  which  was  motivated  by  the  transportation  industry. 

Next  this  thesis  discussed  the  theory  of  interest  to  this  study  and  the  governing 
equations.  The  equations  of  motion  2.2. 1-4  were  stated  along  with  the  boundary  and 
continuity  condition  to  satisfy.  The  coordinate  transformation  was  introduced  and  a 
detailed  description  of  the  dimension  analysis  was  given.  The  equations  of  motion  and 
the  boundary  and  continuity  conditions  were  transformed  and  scaled.  Then  the  steady 
state  system  was  introduced  and  given  2.8. 1-4. 

Chapter  3  described  the  methods  used  to  solve  the  system  and  presented  the 
solutions  to  the  rail  without  coating  problem  as  well  as  the  coated  rail  problem.  The 
Fourier  transfonn  was  introduced  for  the  rail  only  problem  and  an  analytic  solution  was 
found  in  the  transform  domain  by  solving  a  homogeneous  linear  system  of  ordinary 
differential  equations.  The  Fourier  transform  was  inverted  with  respect  to  the  stresses 
rather  than  the  displacements.  The  solutions  in  the  spatial  domain  for  the  no  coating 
problem  are  given  in  3.2.28-30.  Then  a  similar  method  was  used  to  solve  the  coated  rail 
problem  in  the  transfonn  domain.  Since  the  solutions  were  unsuitable  for  inverting  the 
transfonn,  a  finite  difference  method  was  introduced  yielding  the  difference  model  for 
the  equations  of  motion  3.4. 1  and  the  boundary  conditions  3.4.3,  6.  The  chapter  ended 
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with  a  point-wise  comparison  of  the  numerical  method  with  the  analytic  solution  in  the 
transfonn  domain  to  validate  the  code. 

Chapter  4  began  by  illustrating  the  behavior  of  the  material  in  the  rail  as  stress 
propagated  through  the  thickness  of  the  steel  with  respect  to  a  point  source,  a  unifonn 
source,  and  a  parabolic  source.  Next,  the  parabolic  loading  function  was  used  to  examine 
the  coating  problem.  Then  the  effects  of  various  coating  thicknesses  and  shear  moduli 
were  presented.  Finally,  the  effect  of  layering  coatings  was  also  presented. 

5.2  Significant  Findings 

The  focus  of  this  study  was  to  understand  how  coating  thickness  and  shear 
modulus  affect  the  strength  and  distribution  of  stress  in  a  half-space.  This  was  motivated 
by  a  desire  to  mitigate  the  gouging  phenomenon  due  to  high  speed  impact  loading.  The 
coating  thickness  and  shear  modulus  were  limited  by  the  realistic  properties  of  the 
HHSTT.  From  Figures  16-18,  it  was  deduced  that  the  thickness  of  the  coating  in  the 
ranges  of  interest  had  little  to  no  affect  on  the  strength  or  the  distribution  of  the  load  to 
the  rail.  The  lack  of  influence  of  this  property  is  because  the  coating  is  likely  too  thin  to 
see  a  significant  effect  on  the  stress 

Even  though  the  coating  thickness  had  little  effect,  the  shear  modulus  of  the 
coating  did  have  a  significant  effect  on  the  direct  stress,  <j::  ,  in  the  direction  of  motion. 

In  Figures  19  and  20,  as  the  dimensionless  parameter,  G  ,  increases  (shear  modulus  in 
the  coating  decreases),  the  strength  of  the  stress  is  diminished.  This  implies  coatings  that 
are  less  stiff  have  a  significant  effect  on  the  failure  of  the  rail.  Even  though  the  gouging 
phenomenon  was  not  present  in  the  formulation  of  the  problem,  this  is  a  significant  factor 
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to  consider  when  choosing  coatings  to  mitigate  gouging.  The  shear  modulus  does  not 
have  a  significant  effect  on  the  distribution  of  the  load  over  the  surface  of  the  rail. 

It  was  also  shown  in  Figures  22-24  that  having  a  variable  shear  modulus  through 
the  thickness  of  the  coating  has  no  greater  effect  on  the  strength  or  distribution  of  stress 
than  having  a  uniform  coating  with  a  low  shear  modulus.  This  implies  that  layering 
coatings  of  different  material  properties  will  not  have  a  significant  effect  on  stress 
distribution  in  the  rail. 

In  summary,  the  strength  of  the  force  applied  to  the  surface  of  the  rail  can  be 
diminished  by  choosing  a  coating  with  a  low  shear  modulus.  The  distribution  of  the 
stress  over  the  surface  of  the  rail  will  be  virtually  unaffected  by  the  coating  due  to  its 
thickness. 

5.3  Suggestions  for  Future  Study 

This  research  investigated  the  problem  as  if  the  slipper  had  been  traveling  for  an 
infinite  amount  of  time.  It  also  assumed  the  displacements  in  the  rail  and  coating  were 
not  changing  in  the  vicinity  of  the  shoe  with  respect  to  time.  It  would  be  useful  to  see  if 
the  non-steady  state  problem  behavior  was  drastically  different  from  the  steady  state 
solution. 

Even  though  the  numerical  results  can  determine  certain  behaviors,  it  would  be  of 
interest  to  have  an  analytic  representation  of  the  coated  problem  to  understand 
mathematically  why  these  behaviors  exist.  Even  though  the  solutions  in  the  Fourier 
domain  are  unsuitable  for  inverting  the  transform,  approximation  methods  can  be  used  to 
get  an  approximate  analytic  solution.  Since  the  parameter  5  is  so  small,  it  may  be 
interesting  to  perfonn  a  perturbation  and  expand  the  system  3.3.9  in  tenns  of  8  .  It  is 
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expected  that  the  leading  order  behavior  of  the  analytic  solution  to  the  coated  rail 
problem  is  the  same  or  very  similar  to  the  solution  to  the  no  coating  problem. 
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Appendix  A:  Q  Matrix  Derivation  for  Rail  Only  Problem 


When  finding  the  coefficients  for  the  general  solution  to  the  rail  only  problem,  the 
derivations  became  cumbersome  and  the  matrix  Q  was  introduced.  From  equation 


3.2.22, 


0 

q2 


a\p\ 

0 


0 

a\p2 


iar2 


iaS2 
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SF 

-K 


Now,  substituting  3.2.21  in  for  %  produces 
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-ap-i  sgn(-a) 

SF~ 
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Since  \a\ sgn(-a)  =  -a  and  - a  sgn(-a)  =  \a\ ,  divide  the  first  row  through  by  \a\  and  the 
second  row  through  by  -ia .  Also  divide  both  sides  by  G  and  get 


QP  = 


SF 
\a\G 
-iK 
_  aG  _ 


where  the  matrix  Q  is  defined  as 


0  = 
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Appendix  B:  Convolutions  with  Loading  Functions  for  Rail  Only  Solution 


In  3.2.28-30  the  solutions  for  stress  are  written  in  terms  of  convolutions  with  the 
applied  loading  functions.  In  chapter  4,  these  solutions  are  used  to  produce  contour  plots 
to  determine  the  behavior  of  stress  in  the  rail.  In  order  to  use  these  solutions  the  contours 
must  be  performed.  Notice  in  3.2.28-30,  there  are  really  only  two  term  types  to  convolve 
with  the  loading  functions  because  everything  else  is  constant.  They  are 


W,7) 


pW+C 


2 


7 


pW+C 


2 


From  chapter  4,  the  three  loadings  functions  used  are  the  Dirac  delta  function,  the 
Heavyside  function,  and  a  parabolic  function.  First,  notice  that  the  convolution  with  the 
delta  function  is  trivial,  that  is  any  function  convolved  with  the  delta  function  is  just  that 
function.  The  solutions  in  3.2.28-30  are  effectively  the  convolution  with  the  delta 
function. 

Next,  the  uniform  distribution  of  force  applied  to  the  surface  is  accomplished 
through  a  constant  times  the  Heavyside  function,  H(l  -|(T|)  •  In  this  research,  the 


constant  load  applied  was  one  giving  the  convolutions  as 


/*z,  =  \h(  i-|<r|) 
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dx 
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Now  the  Heavyside  function  is  zero  outside  the  interval  -1  <  C,  <  1  so  the  integrals 


reduce  to 


/*Z,  =  J 


C~X  ,  1 

— z^— - pix  =  —  In 

i  (C-xy+(ptrj)  2 


(C  +  l  f  +(p,Tjf 


(c-i  y+(p,uy 


V 


rdx  =  ■ 


1 

tan  1 

U+i'. 

-tan  1 

U-iX 

Pi 

{  Pi7!  J 

{  Pi7!  )_ 

f*Ni=  [ —  , 

i(£-x)  +{Pirl) 

Similarly,  convolving  the  parabolic  load,  4.2.1,  with  the  solution  gives 
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Appendix  C:  Derivation  of  Boundary  Conditions  for  Coating  Problem 


When  finding  the  coefficients  for  the  general  solution  to  the  coating  problem 
system  3.3.9  was  introduced.  The  matrix  is  comprised  of  sub  blocks,  R^ ,  R2,  L  ,  found 
by  using  the  loading  condition  3.3.2  and  the  interface  conditions  3.3.3  and  3.3.4.  Rx  and 
R2  are  derived  by  applying  3.3.3  and  3.3.4  to  the  general  solutions,  3.2.17,  3.2.18,  3.3.6, 
and  3.3.7.  Using  the  interface  condition  3.3.3  at  i)  -  0  produces 

1  1  1  r  1  1  i 

r(2)  (2)  (2)  (2)  A~  .(2)  .(2)  P 

*/V2  *^3  •a'4  J  [_*^i  •a'2 

where  x  symbolizes  the  eigenvectors  associated  with  the  rail  (region  II).  Using  the 
interface  condition  3.3.4  at  tj  —  0  produces 

\  +  iaS2x[2')  X2+iaS1x{2'1  +iaS2x\2'1  A4  +  iaS2x 

iat\  +  ^/l/l1x^2,  iat\  +q^X2x'2'1  iar3  +  qlA3x32'1  iai\+qxA4x 

a\p2 

_<hpV 

s2  J 

Combining  the  two  interface  conditions  together  produces  the  system 

[*>  u[^]=M 

where 
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Now  L  is  derived  by  applying  3.3.2,  at  77  =  1 ,  to  the  general  solutions  to  produce 


[a^  +  iaS2x\2'  j  e*2  {^A2  +  iaS2x  2^  j  eA 3  ^A3  +  iaS2x ^  +  z«<5',x42)  j 

^ iarx  +  qlAlx[ j  +  qtA2x^  j  e ^  |z«f]  +  qtA3x^  j  eA 4  ^zaz-  +  q{A4x^^ 

SF 

-K 


So  let  L  be  the  2x4  matrix  on  the  left  hand  side  of  the  equation  and  combine  with  the 
results  from  above  to  produce  3.3.9: 
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A 

~  0~ 

L 

®2x2_ 

F 
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Appendix  D:  Derivation  of  Boundary  Condition  Matrices  for  the  Numerical  Method 


Equation  3.4.5  gave  the  fonn  of  the  boundary  condition  at  77  =  0  where  H 
captures  the  influence  of  the  rail  on  the  coating.  Recall,  H  is  found  by  substituting 
3.4.4  for  p  in  the  left  hand  side  of  3.2.22  and  substituting  into  3.3.4.  This  produces  the 
following  results: 
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which  with  some  simplification  yields 
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0  \0c\p2 

X 

q2'_ 

Next,  the  difference  equation  was  given  for  the  boundary  condition  which 
produced  a  ghost  point.  This  boundary  condition  was  then  substituted  into  3.4. 1  to  give 
3.4.6  where  the  matrix  G  is 
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